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Abstract 

We study the discrete Schrodinger operator H in Z d with the surface potential of the 
form V{x) = gS(xi) tan 7r(a • X2 + w), where for x G Z d we write x = (x\,X2), x\ G 
Z dl , X2 G Z^ 2 , a G M^ 2 , u; G [0,1). We first consider the case where the components 
of the vector a are rationally independent, i.e. the case of the quasi periodic potential. 
We prove that the spectrum of H on the interval [—d, d] (coinciding with the spectrum 
of the discrete Laplacian) is absolutely continuous. Then we show that generalized eigen- 
functions corresponding to this interval have the form of volume (bulk) waves, which are 
oscillating and non decreasing (or slow decreasing) in all variables. They are the sum 
of the incident plane wave and of an infinite number of reflected or transmitted plane 
waves scattered by the "plane" Z d2 . These eigenfunctions are orthogonal, complete and 
verify a natural analogue of the Lippmann-Schwinger equation. We also discuss the case 
of rational vectors a for d\ = d<i = 1, i.e. a periodic surface potential. In this case we 
show that the spectrum is absolutely continuous and besides volume (Bloch) waves there 
are also surface waves, whose amplitude decays exponentially as \x\\ — ► oo. The part of 
the spectrum corresponding to the surface states consists of a finite number of bands. For 
large q the bands outside of [—d, d] are exponentially small in q, and converge in a natural 
sense to the pure point spectrum, that was found in E9] in the case of the Diophantine 
a's. 
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1 Introduction 



The progress of the last decades in spectral theory of differential and finite difference operators 
with random ergodic and almost periodic coefficients in the whole space makes natural the 
study of operators with same type of coefficients supported on a subspace only. Being of evident 
interest from the point of view of wave physics, they provide a class of operators "intermediate" 
between operators whose coefficients decay in all coordinates (scattering theory) and operators, 
having coefficients of the same order of magnitude in all coordinates. We mention recent papers 
HH, ||, fl9fl , || - [f[4l| ) Hl5| , IfTgfl , devoted to the study of the spectral and related properties 



for operators of such a kind. These operators are either defined on the half-space by random, 
almost periodic or periodic boundary conditions or have the same type of coefficients supported 
on certain subspaces of IR d or Z d . 

As in |19| we consider here the discrete Schrodinger operator 

H = H + V (1.1) 

acting on l 2 (Z d ), where 

(# *)(x) = -1/2 *(y)> C 1 - 2 ) 

|x-y[=l 

is the discrete Laplacian, 

V(x)=6{x 1 )v(x 2 ), x=(xi,x 2 ), ii G x 2 eZ d2 , d 1 + d 2 = d, (1.3) 



with 



v(x 2 ) = gtan7r(a • x 2 + to) 



;i-4) 



is the multiplication operator ("surface" potential), whose support is the subspace Z d2 of the 
space Z d , and 



d u d 2 G N, g > 0, a G R d \ u G [0, 1) 



;i.5) 



are the parameters determining the potential. 

It was shown in Jl9| that for any g ^ 0, u G [0, 1), and for a G M. d2 , satisfying the Diophantine 
condition i.e. there exists e > such that 



\a-x 2 -m\> const /\x 2 \ d2+£ , Wx 2 G Z \ {0}, Vm G Z, (1.6) 

the spectrum of H = H + V, lying outside the spectrum [—d, d] of the discrete Laplacian (|1.2|) , 
is pure point, dense, of multiplicity one, and the respective eigenfunctions decay exponentially 
at infinity. 

The "volume" version of this operator, corresponding to the case di = 0, has been studied 
earlier in . The operator has a complete system of exponentially decaying eigenfunctions, 

corresponding to the pure point dense spectrum of multiplicity one occupying the whole real 
axis. This spectral structure is caused by strong and irregular fluctuations of the quasi periodic 
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potential (jl.3|) . It is the extreme case of the strong localization regime, which in general appears 
either if, for a fixed energy, the amplitude of the potential (random or almost periodic) is large 
enough or, if for a fixed potential, the energy is close enough to the spectrum edges (see [^2 
for related results and references). 

In the case di = 1 the support of the potential is the hyperplane Z d-1 of the space Z d . This is 
why it is natural to call the respective operator ( |1 . 1| ) - ( |1.5| ) the surface Maryland model. This 



operator is closely related to the boundary value problem (|3.5|) , considered in [15], [II], [14| . We 
may also call the operator ( |1 . 1| ) - ( |1.5| ), for d\ > 2, the subspace Maryland model. 

These models can be analyzed in great detail, thereby providing examples of spectral types 
which are only partly known for general random or almost periodic function v in (O). All 



these versions of the Maryland model have an absolutely continuous component of the spectrum. 
This component was first indicated in [15], and then was studied in |Tl| in the context of the 



boundary value problem defined by ( |1.4j ) and by formula ( |3.5| ) below. It was proven that if the 
components of the vector a G M d2 are rationally independent, this part of the spectrum of H is 
purely absolutely continuous and also that the properly defined wave operators corresponding 



to this part exist and are complete. Besides, it was proven in |L4) that the surface states (see 
]T5j, [TJJ for definitions) are absent. 



In this paper we develop several general ideas and results of the theory by considering the 
explicitly soluble model, defined by formulas ( |L1~1) - (|1.5|) . We begin by showing that the 
Green function of the model can be written in a rather convenient form (Section 2). By 
using this form we study first the quasi-periodic case of rationally independent components 
of the vector a in (|1.3|) (Section 3). We prove that the spectrum of the operator is purely 
absolutely continuous on the interval [— d, d] (on the spectrum of discrete Laplacian) and that 
the wave operators, corresponding to this part of the spectrum exist (these facts were proved 
in ]llj, Q by other methods). Then we find an explicit form of the generalized eigenfunctions 
(polynomially bounded solutions of the respective equation), corresponding to this part of the 
spectrum. These eigenfunctions possess properties, similar to those of the Sommerfeld solutions 
of scattering theory. Along the x<i direction, they behave like Bloch-Floquet solutions. They 
are orthogonal and complete on the interval [—d, d] of the spectrum. As they do not decay 
in the longitudinal coordinates x\ we call them volume states. We consider also the case of 
rationally dependent components of the vector a in ( |1.4| ), where the respective surface potential 
is periodic in x 2 > restricting ourselves to the technically simplest case of d\ = = 1 (Section 4). 
In this case the whole spectrum is absolutely continuous. It consists of the interval [—d, d] as in 
the quasiperiodic case, and of a certain number of intervals, some of them possibly intersecting 
[— d, d]. To the interval [— d, d] correspond generalized eigenfunctions which do not decay in 
the longitudinal coordinates x\. Instead the generalized eigenfunctions corresponding to the 
other intervals decay exponentially in Xi, being of the Bloch-Floquet form in the longitudinal 
coordinate x 2 . Such a type of surface states (see Definition |37l] below) were first found by 
Rayleigh in the problem of oscillation of an homogeneous elastic half-space (see e.g. |]2l|), and 
since then were found and studied in a number of problems, described by differential and finite 
difference equations whose coefficients are strongly varying in coordinates x\ (see e.g. 0] for 
a list of references on respective physics results and applications). All these results concerned 
the problems where coefficients were constant in the X2 coordinates. We analyze also the case 
where a n = p n /q n approaches an irrational aasn-> oo, and we show that there exists a certain 
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continuity of the spectrum in this asymptotic regime. In particular, the width of surface bands, 
lying outside of [—d, d] is exponentially small in q n as n — > oo, and the bands approach the 
dense set of eigenvalues, found in |L9| . 



2 Generalities 

Recall that we are studying the self-adjoint operator H, acting in Z 2 (Z d ) and defined in ( |1 . 1| ) 
- (|1.5|). The operator is selfadjoint as the sum of the multiplication selfadjoint operator V of 



flOp, and of the bounded selfadjoint operator H of (|1.2|) . We will use an analogue of the 
Cayley transform introduced in for the "volume" potential (d± = 0) and in Jl9| for the 
"surface" case (di = 1), in both cases to study the pure point spectrum for the Diophantine 
a's (see (Ofi). 



To put the subsequent simple argument in a more general context, we rewrite the potential 
(O) as 



V{x)=v(x 2 ) X s(x), (2-1) 

where xs is the indicator of the subspace S = Z d2 and we assume that g > (the case g < 
can be treated analogously). We define the orthogonal projection P of l 2 (Z d ) : 

(P$)(x) = X5 (x)$((0,x 2 )), (2.2) 

and we write the potential ( |2.1| ) in the form 

V = PvP, (2.3) 

Here and in the following we use lower cases to denote operators acting on l 2 (S) defined by the 
restriction on Pl 2 (Z d ) of the corresponding operator. 

We use as a starting point the well known formulas for the resolvent G(z) = (H — z)~ l of a 
selfadjoint operator H = H + V: 

G{z) = G (z) - G (z)T(z)G Q {z), %z + (2.4) 

with 

G (z) = (H - z)-\ T(z) = V- T(z)G (z)V. (2.5) 
It follows from (|2.3|) and from (|2.5| ) that the operator T(z) has the form: 

T{z) = Pt{z)P, (2.6) 
where the operator t(z), acting on l 2 (S), satisfies the equation 

t(z)=v-t(z)>y Q (z)v, (2.7) 
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in which 70(2) is defined from the restriction of Gq{z) to the subspace P/ 2 (Z d ). The formal 
solution of the equation is 

t(z) = v(l + ^vy 1 = (V 1 + i (z))-\ (2.8) 

Let u be the unitary operator in l 2 (S) defined by the relation: 

(uiP)(x 2 ) = e- 2i ™- X2 iP(x 2 ), x 2 E S. (2.9) 

Then, by using the Euler formula for the function x i— > tanx and the notations above, we can 
write the potential Q1.4J) as 

q 1 — au , 

v = ^-——, 2.10 
1 1 + au 

where 

a = e - 2inuJ . (2.11) 



Formulas (|1.3|) - Q2.11 ) motivate the following abstract statement. 



Lemma 2.1. Let H be a selfadjoint operator, acting on l 2 (7L d ), and having the form H = 
H + V , where H is a selfadjoint operator and V is given by formulas l \2.$ )and( \2.1(\ ) in which 
S is any subset of Z d and \a\ < 1. Define the following operators in l 2 (S) 

b(z) = (g l0 (z)-i)(g l0 (z)+tr 1 , (2.12) 

assuming that b(z) is bounded. If the operator g'yo(z) + i is invertible and if 

\\b(z)\\<l, (2.13) 

then the operator t(z) , defined in ( \2.6\ ) and in can be represented in the form: 

t(z) = g(l - au)(l - ab{z)u)-\g l0 (z) + z)~\ (2.14) 

or in the form 

q-l 

t(z) = gig^+i)- 1 ^ - 2iau^2{ab(z)u) l (l - {ab{z)uY)- x {g lo {z) + i)' 1 } , (2.15) 

1=0 



where a is defined in ( 2.11 ), and q > 1 is an integer. 



Proof. Note that the conditions < 1 and \o~\ < 1 allow us to define the operator 

(1 — cr6(z)-u) _1 by the Neumann-Liouville series. Consider first the case, where the modulus of 
the complex number a in Q2.10D is strictly less than 1. In this case the operator (1 + au)' 1 is 
well defined and we obtain from ( p.!0| ), and from ( |2.12| ): 

l + 7o^ = + o~u) + #7o(l - au)] (i(l + au)) 1 

= (#70 + (91o - i){gio + i)~ x ov) + o-u))' 1 , 



or 1 +7of = (<77o + «)(l — b(z)au) + o~u)) , where the operators 70(2), and b(z) are defined 
in (|2.12| ). Formulas (|2.7|) , ( p. ICQ , and the hypotheses of the lemma lead to (|2.14|) for \a\ < 1. 
According to inequality fl2.13| ) the Neumann-Liouville series for (1 — b(z)au)~ l converges for 
I a I = 1, and since the operator (1 + au)~ l is not present in formula ( |2.14| ), we can make the 
limit I a I — > 1 in the formula, proved for |cr| < 1, and obtain representation ( 2.14j ) in the case 
lo-l = 1. 



Proposition 2.1. Let H be the selfadjoint operator defined in Lemma \2J\ and G(z) = (H — 
z) -1 , Qz > be its resolvent. Assume that z is such that the conditions of the Lemma pO| hold. 
Then G(z) can be represented as follows: 



G(z) = G (z) - gG {z)P(g lo (z) + i)- l PG a {z) + 2igG (z)P(g l0 (z) + i) 



q-l 



o-uJ2(°Kz)u) l (l - {ab{z)uyy\g lQ {z) + iy 1 PG (z), 



(2.16) 



1=0 



where q > 1 is an integer, u is defined in and the operators 70 {z), b(z) are defined in 



Proof. The proposition follows easily from (|2.4j) , and from Lemma ^2A . 



Remarks. 1). In formula ( |2.10| ) the unitary operator au can be viewed as the Cayley transform 
of v (see [|l| for the definition of the Cayley transform). Likewise, the contraction operator b(z) 
can be viewed as the Cayley transform of the dissipative operator ijo(z) (Jfti'y > 0). Hence, we 
can say that the passage from the operators v" 1 and 70(2) in (|2.8| ) to their Cayley transforms 
au and b(z) in the case of the potential ( |1.3| ) - ( |1.4| ) leads to formulas ( |2.14D - ( |2.16| ). This 
will allow us to study the absolutely continuous spectrum of the operator H for any d\ > 0, 
as it was done in papers (22| and [19| for the pure point spectrum, despite that the subsequent 
techniques to study the resolvent ( [2.16| ) are different in these two cases. 

2). Integrate formula ( p,16|) with respect to uj G [0, 1) and denote this operation by (•••). We 
obtain: 

(G(z)) = G (z)-gGo(z)P(g lo (z)+t)- 1 PGo(z). 



In view of the general formula ( p.8|) , valid for any surface potential v, we can interpret the 
equality (t(z)) = gig^oi 2 ) = (—(w)" 1 +7o( 2 0)~ 1 as the fact that (G(z)) is the resolvent 

of the Schrodinger operator whose surface potential is the complex constant V(x) = —igxs{ x )- 
This fact plays an important role in the interpretation of results of analysis of the point spectrum 
of H outside [— d, d] in Similar fact is known also in the case of the volume potential ( |2.3| ), 
i.e. for the case S = Z d §. 

Now we are going to show that the above proposition is applicable to the operator defined by 
(|1.1| ) - (|1.5|) where S is chosen as Z d2 . To check the conditions of the lemma and the proposition 
we will use the Fourier transformation which we define as follows: 

$(k) = e~ 2mx ' k $(x), k e T u , $(x) = f dke 2i " x - k l>(k), x G Z v , (2.17) 
where T v = [0, l) u is the i/-dimensional unit torus. 
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By using the Fourier transformation we can write the following representation of the Green 
function Gq (x — y; z) of the z/-dimensional Laplacian (operator (|1.2|) for d = v): 

<%\ic-y ] z)= [ dk 6 ^ ( * Sz^O, (2.18) 

where 

E u (k) = -^cos27rA; i , (A*, fc„) = k e T v . (2.19) 

i=l 

These formulas imply that the operator 70(2) of (|2.12| ) has the following matrix in / 2 (Z (i2 ): 

lo(x 2 -y 2 ;z) = G { d \(0,x 2 ) - (0,y 2 );z), (2.20) 
i.e. 70(2) is a convolution operator in Z 2 (Z d2 ). In view of (|2.18| ) we have: 

(2 - 21) 



T d 2 



or 



7o(x 2 ;^)= f dk 2 e 2i ^G^(0,z-E d2 (k 2 )). (2.22) 

We shall denote 

%(k 2 ; z) : = G ( ^\o, z - ^(fe)), (2.23) 

i.e. %(k 2 ;z) is the symbol, representing the operator 70(2) in L 2 (T d ' 2 ) as a multiplication 
operator. These formulas allow us to show that the hypotheses of Lemma |2.1| and Proposition 



2.1 are valid for any z, ^sz > (see Lemma |5T2j ) . Besides, we have 

Lemma 2.2. Letb(z) andu be the operators, defined by ) and fljj.fy ). Then for any integer 
m > I, 



m—l 



((b(z)u) m ^){k 2 ) = a m (Y[b{k 2 + la;z))0{k 2 + ma), k 2 eT d \ (2.24) 



1=0 

where ip denotes the Fourier transform of <p e l 2 (Z d2 ) and 

kk,-,z) = ^'\- i (2.25) 

where ^ (k 2 ,z) is defined in (\2.23j ). 

Proof. It follows from ( |2.9|) that the operator u is the shift by a in the space L 2 {T d2 )\ 

(wp)(k 2 ) = 0{k 2 + a). (2.26) 

From this and the fact that b(z) of (|2.12f) is the multiplication by the function b(k 2 ; z) of (|2.25| ) 
in the space L 2 (T d2 ) prove the lemma. 

We will obtain now a representation of the Green function of H which will be central in the 
subsequent spectral analysis of the absolutely continuous spectrum of the operator. 
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Theorem 2.1. Let H be the operator, defined by ( |i.i| j-( [7T5| . Then the Green function of H 
(the matrix in l 2 (Z d ) of its resolvent G(z) — (H — z)~ x ) can be written in the form: 

oo „ 

G(x,y;z) = G^x-y;^^ dk 2 e 2 ^ x ^t m {k 2] z) 

x Gt\x x] z - E d2 (k 2 ))G { dl \ yi ; z - E d2 {k 2 + ma))e- 2 ™ m ^\ (2.27) 

where 



-1, m = 

t m (k 2 ;z) = - - - - ■ , { 2ia(g%(k 2 + a;z) m = 1 

2ia m+ \g%(k 2 + ma-z)+iy l ]]™- 1 b(k 2 + la-z), m>2, 

(2.28) 



g%(k 2 ;z) +i 



Gq (x\;z) is the Green function &2.1Q ) of the d\- dimensional Laplacian, E d2 (k 2 ) is defined 
( \2.1Q) for v = d 2 , and 70(^2; z), b(k 2 ; z) are defined respectively in ^2. 2S(j and ( \2.23j) . 

Besides, the (generalized) kernel of the operator T(z) of ( \2.4\ ) and of Lemma \2.1\ has the fol- 
lowing form in L 2 {T d ): 

00 

T(k, p; z) = t m {h; z)8(k 2 + ma- p 2 ), (2.29) 

m=0 

where t m [k 2 \ z) is defined in ( \2.2fy ). In particular, the kernel is independent of the components 
k±,p\ G T dl of its arguments k, p e T d . 

Remark. Formulas (p.27|) and (|2.29| ) have to be compared with the formulas for respective 
quantities for point potential: V(x) = t>5(x), (d 2 = 0) and for the constant surface potential: 
V(x) = v5(xi), v = const. In the first case we have: 

G(x, y; z) = (x - y; z) ^— G?'(x; z)G^(y; z), (2.30) 

1 + vG ( d) (0;z) 

and 

T(k,p;z) = " (2.31) 

l+vG { o } (0;z) 

while in the second case: 

p 2iwk 2 -(x 2 -y 2 ) 



III 



G(x, y; z) = G^ (x-y;z)-v / dk 2 - 

JT d 2 



and 



h<h ~l + vG^\0;z-E d2 (k 2 )) 
x Gi? 1 \x 1 -z-E d2 (k 2 ))G {dl \y 1 -z-E d2 (k 2 )), (2.32) 



t/i \ vS(k 2 -p 2 ) 

T(k, p; z) = 7^-. . (2.33) 

l + vG {dl) (0;z-E d2 (k 2 )) 
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In particular the term, corresponding to m = in fl2.27j ), coincides with the second of ( 2.32|) 
in which v is replaced by ig. 

Proof of Theorem (\2.1\) . According to (|2.13|) , ||6(^)|| < 1 if Qz ^ 0. Hence we can write the 
operator (1 — abu)^ 1 in ( j2.14| ) for q = 1 as the Neumann-Liouville series in powers of abu. 
Applying lemma to each term of the series , we get ( p. 27 ) after simple algebra. Formula 
( |2.29| ) follows from ( |2.4j ) and fl2.27| ). Theorem |2.1| is then proved. 

Remark. Formulas (|2.27|) and Q2.29| ) are the basic tools of spectral and scattering analysis of 
the operator ( |1.1| ) presented in this paper. An advantage of these formulas is that they are 
valid for all values of the spectral parameter z = E + ie, up to the real values z = E ± zO, 
for \E\ < d, in the case of a's with rationally independent components (quasi-periodic in x 2 
potential V(x)) and they are valid for all E G R in the case of a's with rational components 
(periodic in x 2 potential V(x)). 

One more general fact, concerning the operator H and necessary in the sequel, is given by 

Theorem 2.2. Let H = H + V be the operator defined by ( [Op , ( jOU and ( pHp . Then its 
spectrum cr(H) contains the interval [—d, d] = o~(H ) for all g El, a G M. d2 and uj G [0, 1]. 

Proof. We will apply the H. Weyl criterion, according to which £eR belongs to the spectrum 
of a self-adjoint operator H if and only if there exists a sequence {\I/ n } ne N of vectors of respective 
Hilbert space such that \\^? n \\ = 1, and that lim^oo \\(H — E)ty n \\ =0. 

Denote by l r the indicator of the ball {x G Z d : |x| < r} and set for all k G T d , 

*„(x) = l n (x)(l - 5(x 1 ))e 2OTk - x /iV„; N 2 n = ^ |l„(x)(l - 5{ Xl ))\ 2 = 0(n d ), n - oo. 

It is easy to find that 



(^*n)(x) 



£ d (k)^ n (x), |x| <n-2,\ Xl \ >2; 
A n (x), n - 2 < |x| < n + 2; 

6n(x), < 1; 



0, |x| > n + 3, 

where \\A n \\ = 0(n _1//2 ), ||6 n || = 0(n~ dl / 2 ) as n — >• oo. This proves the theorem. 



3 Absolute Continuous Spectrum in the Almost Periodic 

Case 

In this section we assume that the vector a G M. d ' 2 from (|1.3| ) has rationally independent 
components, i.e. that the relation a\r% + ... + ad 2 Td 2 — with rational coefficients rx,...,rd 2 
implies that all these coefficients are equal to zero. 

Theorem 3.1. Let H = H + V be the self-adjoint operator defined by - ( \L.dJ in which 
the vector a G M. d2 has rationally independent components. Then H has purely absolutely 
continuous spectrum on the interval (—d,d). 



9 



Proof. According to the general principles (see e.g. pi[), it suffices to prove that for any 
vector $ G Z 2 (Z d ) of a dense set the limit Q(G(E + i0)$, $) exists and is bounded for all 
E G (—d,d). Restricting ourselves to the vectors concentrated at a point x G Z d , i.e. to the 
vectors <5 X = {5(x — y)} ye z d , we have to prove that for any x G Z d the limit 3G(x, x; E + z'O) 
exists and is bounded for all E G (—0?, d). We shall prove more, namely that G(x, y; E + z'O) 
exists and is bounded for all E G (— d, d) and all x, y G Z d . In view of Theorem we have 
to prove that the series of (|2.27|) converges not only for > but also for Qte = 0. 

Since the vector a has rationally independent components, we have uniformly in fc 2 and 
for any 7 > (see e.g. [[|): 



lim |{leZ: k 2 + lae K J (E), 1<1< mjm" 1 = |AT 7 (S)| 



(3-1) 



where 



K,{E) = {k 2 G T d ' 2 : £ - £, 2 (£; 2 ) G + 7, <*i - 7]}, 



(3.2) 



and \K y (E)\ denotes the Lebesgue measure of the set K^(E) C T d2 . It is easy to check that 
for any \E\ < d there exists 7 > such that K y (E) is an open set of T d2 . According to Lemma 
p73| , in this case there exists 5 > such that |6(/c2, E + i0)\ < 1 — 5, \/k 2 G K y (E), and according 
to (|3.1|), there exists m > such that 



m , 



G Z : k 2 + lae Kj(E), 1 < I < m} > —\K^(E)\ 
for all m > m . Hence we have the following bound for the product in the r.h.s. of ( p. 27] ): 



m—l 



[ J b(k 2 + la; E + i0) 



1=0 



<(l-5) 



m\K^E)\/2 



m > mo, 



(3.3) 



and the series in the r.h.s. of fl2.27|) converges uniformly in k 2 G T d ' 2 . Besides, by using bound 
(|3.3|) and Lemma |5.5| , it can be shown that for \E\ < d — 7, 7 > 0, the series is bounded in k 2 
and , hence we can integrate the series with respect to k 2 . Theorem is proved. 

Remarks. 1). Another form to express (|3.1|) - (|3.3|) is to write the relation: 



lim 

m— >oo 



m—l 



\\ b(k 2 + la;E + z0) 



i=0 



l/m 



exp 



dg 2 log |6(g 2 ; ^ + *( 



T rf 2 



(3.4) 



valid uniformly in k 2 G T d2 (see 0) and showing that if \E\ < d — 7, 7 > 0, then the integral in 
the r.h.s. is negative, thus the product in the l.h.s. is exponentially decaying in m as m —>■ 00. 

2). Theorem [3.i| reveals a fairly simple mathematical mechanism responsible for the absolutely 
continuous spectrum for the "subspace" potential (|1.3Q - (|L4 ) with di > 1 (recall that in the 
"volume" case d\ = 0, d 2 = d, the absolutely continuous spectrum is absent, moreover if a is 
Diophantine then the spectrum is pure point ||). The mechanism is the positiveness of the 
imaginary part of 70(^2; E + iO) = G Q dl \o,E + iO — E(k 2 )) in a certain domain of (E,k 2 ). 
This is most transparent in the "genuine surface" case di = 1, where (0, E + zO) is pure 
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imaginary if \E\ < 1 and is pure real if \E\ > 1, (see formula ( |3.37| ) below). In the latter case 



\b(k 2 ; E + i0))\ = 1 and the series ( p. 27] ) diverges for a dense set of energies (see |1|). This 

leads to the pure point spectrum everywhere outside of the spectrum a(H ) of the Laplacian 

(similarly to the volume case p], where the analogue of 70(^2! E) in ( p.25| ) is real for all 

In the former case | &( A; 2 ; E + ^0))| is strictly less than 1 for any E G (— d, d) on an open set of 

k 2 G T d \ the series is convergent and the spectrum inside of u(H ) = [—d, d] is pure absolutely 

continuous. 

As usual in scattering theory, a fact of primary interest is the existence and completeness of 
wave operators Q± = s • lim^ ±00 e ltH e~ %tH ° S (A) , where £ is the resolution of identity of H , 
and A is an interval of the spectral axis. In the next theorem we prove these properties in our 
case. 



We mention first that in papers || ||, |TT|, the scattering theory was developed for the 
operator Hi, acting in / 2 (Z^_), = {(xi,x 2 ) G 7L d \x\ > 0,x 2 G Z d_1 }, and defined as: 

1 1 )[ ) \ x 2 ) + J2\ X2 - V2 \=i * (0, V2) + v(x 2 )^(0, x 2 ), xi = ^ 

for certain random and almost periodic surface potentials v. The operator can be viewed as 
a boundary value problem for the discrete Laplacian in Z 2 (Z+) with the boundary condition 
\l/(— l,x 2 ) = v (x 2 )^/(0, x 2 ), x 2 G Z d_1 . The "unperturbed" operator if here is the discrete 
Dirichlet Laplacian, corresponding to v = in ( |3.5| ). The operator ifi is closely related to our 
operator H of (|1.1| ) for the surface case d\ — 1, d 2 — d — 1 via standard Green's formulas. 



Theorem 3.2. Under the conditions of the Theorem i| , £/ie wave operators Cl± for the pair 
(H,H ), defined by (fITl|j- fli-fl ), exzsi and are complete for any closed interval A = [a,b] C 



(— cf, d). 

Proof. Existence of wave operators is a rather general fact. It was proved in M for a general 
surface perturbation v in ( |3.5|) . In our case the proof is practically the same. Thus we have to 
prove the completeness. Mimicking the argument of |TT|, , developed for the boundary value 
problem ( |3.5|) , it is easy to reduce the proof of completeness to the proof of the relation: 

sup V \G((x u x 2 ),y;E±ie)\ 2 < 00 (3.6) 

e>0,Ee[a, b ] X2&d2 

for any fixed x x G Z dl , y G Z d and [a, 6] C {-d,d). Our formulas (|2~2"T1 ) - (|2~2"gl ) for the Green 
function of H can be written in the form: 

G((x u x 2 ),y;z) = [ dk 2 e 2 ™ k ^G{(x u k 2 ),y-z), 

JT d 2 

where 

00 

G((x u k 2 ),y;z) = G^'fa - yi, z - E d2 (k 2 )) - ^2t m (k 2 , z)G { dl) (xi, z - E d2 (k 2 )) 

x G^iyi, z - E d2 {k 2 + ma))e 2inV2 - {k2+ma) . (3.7) 



11 



Thus, applying the Parseval equality for the Fourier transform with respect to the variable x 2 , 
we can present the sum in ( |3.6| ) as: 



I ^ 2 |G((x 1; A; 2 ),y;£ + ^)| 2 . 



(3i 



We have shown in the proof of Theorem |3.1| that the series ( |3.7| ) converges uniformly in k 2 G T d2 
for z = E + ie,E G [a, 6] C (— d, d); e > 0. Hence the integral in ( |3.8|) is finite for these values 
of E and e. This proves (^.6[). 

In the next theorem we construct a family of generalized eigenfunctions of H, relating them to 
the Green function of the operator, as in the conventional scattering theory |23|, . 



Theorem 3.3. Let G(x, y; z) be the Green function of the operator H = Hq + V , defined by 
( |OP - (\l-3j) , in which the vector a is rationally independent. Set 



g(x, k; z ) = Y, G ( x > y; ^ k ' y , kGT d 



yez" 



tt,(x,k) = (£ d (k)-z)G(x,k; z) 



and 



(3.9) 



(3.10) 



t d2 = T d2 \ {(0,0,...,0),(7r J 7r,...,7r)}; T d = T dl x T a 
Then, for z = £><f(k) =F limits: 



(3.11) 



#±(x, k) = hrn W*(x, k) | 2=Bd(k)Tie = Hm ±ieG(x, k; £ d (k) T (3.12) 

exist /or a// k G T d , are bounded in x G Z d /or any k G T d , are continuous in k varying in any 
compact set of T d , and have the form 



tt±(x, k) = e 2mk x + ]T t m {k 2 - ma; z)G^ 1 \x 1 ; z - E d , 2 {k 2 - ma)) 



m=0 



2in(k2 —ma)-X2 



=£ d (k) T iO 



(3.13) 



where the coefficients t m (k 2 , z) are defined in 
Moreover: 

(i) the functions \&±(x, k) satisfy the Schrddinger equation in x for any k G T d : 

((^o + V)* ± )(x > k) = S d (k)* ± (x,k); (3.14) 



fiij the functions 1 I r ±(x, k) are tae unique solutions of the equation: 

*±(x, k) = e 2 ^ k x - ^ Gf (x - y; £ d (k) =F i0)^(y)*±(y > k). 



(3.15) 
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for anyk G T d in the class of sequences ^ = { , I'(x)} xgZ d whose restrictions ip = {^(0, x 2 )} X2& d 2 
and the sequences {(1 + ae~ 2l7Ta ' X2 )ifj(x 2 )} X2&Z d 2 are representable as the Fourier transforms of 
measures of bounded variation on T d2 , and the sum of the r.h.s. of ( \3.13i ) is understood as the 
generalized convolution of respective functions and measures; 

(Hi) the families {^±(-, k)} kg fd are orthonormalized, i.e. if for any continuous function $ of 
compact support in T d we set: 



$ ± (x) 



tf±(x,k)$(k)dk, 



then for any two such functions $W and $( 2 ) we have: 



^^(xVl'V (X) 



dk<f>«(k)$( 2 )(k); 



(3.16) 



(3.17) 



iv) the functions \l/± : 7L d x T d — > C are i/ie kernels of the wave operators Q±, whose existence 
and completeness are proved in Theorem \3.3j, i.e. for any $ G l 2 (Z d ) such that the support of 



its Fourier transform $ is a compact set in T we have: 



(n±$)(x) 



tf±(x, k)$(k)dk. 



(3.18) 



Proof. We use again our basic formulas (|2.27 ) - ( |2.28|) for the resolvent of H. Making the 
Fourier transform of ( 2.27|) with respect to y and multiplying the result by Ed(k) — z, we 



present (|3.10|) in the form: 



**(x, k) = e 2mkx + tmifa - ma; z)G^ 1 \x 1 ; z - E d2 {k 2 - ma))e 



2in(k2—ma)-X2 



(3.19) 



m=0 



Each term in this series is continuous in k 2 and z, Qz > and can be extended to real z = 
E + iO,E G [a, b), if the closed interval [a, 6] lies strictly inside (— d, d). According to bound 
( |3.3j ), the series converges uniformly in k 2 G T d2 and a < 9ftz < b,Qz > 0, hence it defines a 
continuous function in this domain. This allows us to perform the limits (p. 12 ) for k 2 G T d2 
and to obtain formula ( |3.13|) . 

Our limitation k G T d , where T d is defined in (|3.11| ) is necessary because for k G T d \ T d and 
for the respective two values of E = ±d we cannot guarantee the validity of bound (|3.3|), thus 



the convergence of the series in formula Q3.13|) . 

Let us prove now property (i) of \&±(x, k). We have obviously: 

#o(x - t)G(t, y; z) - EG(x, y; z) + V(x)G(x, y; z) = zeG(x, y; z) + 5(x - y) 

It— x|=l 



The definition ( p. 10| of ^(x, k) and an easy justification of the interchange of the multiplication 
by V(x) and of the Fourier transformation in y in the third term of l.h.s. lead to the equality: 

£ #o(x - t)tt,(t, k) - £*,(x, k) + V(x)tt,(x, k) = z£(^ z (x, k) + e 2 - k ' x ). 

|t-x|=l 
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Now, in view of relation ( J3.12 ) , the limit of the r.h.s. of the last equality is zero as e — > 0, and 
we get flSrni) . 



Let us prove now assertion (ii) of the theorem, i.e. that \l/±(x, k) satisfy (the Lippmann- 
Schwinger) equation fl3.15| ). We remark first that any solution ^ of ( [3.15|) is uniquely determined 
by its restriction ip(x 2 ) = ^((0, x 2 )) to the subspace Z d2 , and that ip verifies the equation, that 
can be symbolically written as: 



$(x 2 ) = e 2tnk2 ' X2 - (70^)^2), x 2 G Z d2 . (3.20) 

Hence we have to verify that the restriction i[)(x 2 ,\c) of ( |3.13| ) to if 2 satisfies ( |3.20| ). By using 
( |3.19| and ( |3.13|) , we can write the restriction symbolically in the form: 



V = {1 - (<?7o + i)~Vyo[l - 2iau(l - ahi)-\g lG + i)~ 1 ]}e 2 , (3.21) 

where 62(22) — e 2tnk2 ' X2 and we used definition (|2.12|) of 70. The symbols 70, b and u in the 
formula denote now not operators on Z 2 (Z d2 ) or in L 2 (T d2 ), defined in (p32|) and in but just 

operations acting on sequences (functions of x 2 G Z d2 ) and representable as Fourier transforms 
of measures of bounded variation depending on the parameter z = E(k) =FzO, k G TP 2 . In order 
words, they belong to the linear manifold: 

Ac = {f(x 2 ),x 2 G Z d2 : f(x 2 ) = [ e 2i ^ 2 - X2 M k (dp 2 ); VarM k < 00}. (3.22) 

JT d 2 

The operations b and 70 are multiplications of Mk by b(p 2 , z) and by 70(^2, z) with z = ^(k) =F 
iO, and u is the shift by a of the measure. The operation (1 — bu) _1 is defined by the series 
'Y^ =Q {bu) m ', whose terms are given by (|2.24j ), and which converges for all k G TP. By using 
these facts and a simple algebra, we can rewrite ( [3.21D as: 



%l) = i(l + au)(l-abuy\g^ + i)- 1 e 2 . (3.23) 

z=£ d (k)=FiO 



Hence we have for the r.h.s. of (13.201: 



e 2 - 'jovtp = e 2 - 07 O (1 - au){l + au) 1 (1 + au)(l - abu) 1 (#7o + «) ^2 



or 



e 2 -7o«V' = { 1 -fl r 7o(l-<"*)(l-^6") ^To + i) 1 }e 2 , 
meaning that the complex spectral parameter z is replaced by E(k) =p zO. The r.h.s. of the 



relation coincides with ip. To prove this fact we have to repeat the arguments leading to (|3.13| ) 
and ( |3.21| ), but starting from formula ( |2.14| ) for the operator T(z) instead formula ( |2.15| ) . Thus 
we have proved that (|3.13| ) solves ( |3.15|) . 



Let us prove that ( |3.13| ) is the unique solution of ( |3.15| ) in and such that their multiplication 



by (1 + ae 2mh2 ' x ' 2 ) belongs also to C^. Consider the homogeneous equation, corresponding to 

X = lovx (3.24) 
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on the same manifold, and write the equality x — (1 + u )Vi where <p also belongs to ( |3.22| ). 
Then we obtain the following equation for ip: 



(1 +gj )(l - abu)ip = 0, 



where the symbols 70, b, and u are again understood as operations in the class £ k - Applying to 
this relation the operation (1 — abu)~ 1 (g'yo + i)~ 1 ,which is well defined in £ k , we obtain: ip = 0. 

According to the above considerations the second term in the r.h.s. of ( ^.15[ ) is the Fourier trans- 
form of the product of Gq(xi\ E d (k)—E d2 (k 2 )+i0) (the Fourier transform of Gq (x; z)\ 
in x 2 ) and of the measure M k , corresponding to vip: 



\z=E d (\i+iO 



J d 2 



G^\x i; E d (k) - E d2 (p 2 ) + «0)M k (dp 2 ). 



(3.25) 



In view of (I2.10Q , and (|3.23|) we have: 

vip = g(l -au){\ - abu)~ 1 {g^ + i)~ l e 2 



g{gio + i 



-1, 



E d (k)=FiO 
00 

2iY,™{abu) m {g l0 + i)- l e 2 ) 



m=0 



z=E d (k)TiO 



(3.26) 



By using this relation and the notations introduced in Lemma |2.1| and in Theorem [LA 
obtain that the measure corresponding to the second term of the r.h.s. of ( [3.15|) is 



we 



} j t m (k 2 - ma; E d (k) ± i0)5(k 2 - ma - p 2 ) 



m=0 



Combining these formulas we obtain (|3.15|) . 

Let us prove now the orthogonality of \l/±(x, k), corresponding to different k's, i.e. relation 
( |3.17| ). It is clear that it is sufficient to prove ( |3 . 1 7| ) for $W = $( 2 ). The proof is rather 
technical and we outline only its scheme, considering, say 

The first step is the proof of the relation: 



0. 



(3.27) 



where (cf (|3~lp ) 



^ d(k)+i£ (x,k)<l(k)dk, 



(3.28) 



and \l/ z (x, k) is defined in ( |3.10| ), i.e. ^/ E^+iei^, k) = — ieG(x, k; E d (k) +ie). The proof is 
based on formulas ( |3.i9| ), and ( |3.13| ), and on the continuity of G^\x.,E + is) with respect to 



e > 0. It is given in Lemma |3.1| below. 
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The second step is the proof of the relation: 

lim £|<f> £ (x)| 2 = / l*(k)l 2 <fc, ( 3 - 29 ) 



which implies ( fj.!7| ). We will just sketch a proof of this relation. 
Write the resolvent identity for the pair G(z') and G(z"): 

£ G(t, y; z")G(t^;z') = (z' - z")- 1 (G(t, x; z') - G(t, y; z")) . (3.30) 



Replace in the r.h.s. of the identity G by Go — GqTGq (see (|2.4j) ). We obtain after a simple 
algebra: 

G' Gl + {z' - z")- x {G' Q T'G' Q - G'XG'f), (3.31) 

where G' = (Hq — z')' 1 , G' ' = (Hq — z")~ l and T" and T" are the T-operators for the spectral 
parameters z' and z" respectively. Now we make the Fourier transformation with respect to x 
and y, multiplying ( |3.30| ) and ( |3.31|) by e 2l7r p y- 2i7rk x anc i summing the result over x, y 6 Z d . 
The l.h.s. of the obtained relation is (ty z », ^ z ')- As for the r.h.s., it can be written symbolically 

as: 

*(k - p) - T( k , p; (—^ - ^) - T( k , p; ,») (^L_ _ , 

(3.32) 

where T(k, p; z) is the kernel in L 2 {T d ) of the T-operator, whose expression is given in Theorem 
fT[]. Setting in ( P?2D , ^' = E d (p) + ie, and z" = E d (k) + ie, we obtain: 



tf(k-p) - (T(k, p; E d (k) + ie) + T(k, p; E d (p) + ie^ 



x 



E d (p) - E d (k) + ie E d (j>) - E d (k) + lie) ' 



After multiplication by $(k)$(p), where $(k) is a smooth function whose compact support is 
strictly inside T d , and after the subsequent integration with respect to k, p G T d , the second 
term of the last expression tends (weakly) to zero as e — > 0. We use the explicit form ( |2.29| ) 



of the kernel T(k, p; z) to prove that T(k, p; E + ie) is weakly bounded in e > 0, if k, p are 
strictly inside of T d and \E\ < d. After that we are left to prove that the expression in the 
parentheses tends weakly to zero as e — > 0. This proves assertion (iii) of the theorem. 

Let us prove assertion (iv), according to which the solutions \l/±(x, k) are the kernels of the 
wave operators Q±, whose existence and completeness are proved in Theorem |3.2] (see also [14 
for similar results). We will base the proof on the formula (see |[23|| ): 



/oo 
G(E + ie)£ (dE)^, 
-oo 
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where £ is the resolution of identity of the Laplacian H Q of (|1.2|) , and G(z) — (H — z) 1 . In 
the (x, k) representation, usual in the scattering theory, this formula can be written as follows: 



where 



(Ot$)(x) = s - lim (fi e $)(x), (3.33) 

(fi 6 $)(x)=/ ^ d(k)+i£ (x,k)l>(k)dk, (3.34) 
and \l/ 2 (x, k) is defined in (|3.10| ). 

According to general principles, it suffices to prove (|3.33|) for a dense set of vectors of L 2 (T d ). 
We choose a set of functions of the form <E>((&i, /c 2 )) = $1(^1)^2(^2), where $1^ are smooth 
and the support of $1 does not contain the critical points of E^. Denoting the r.h.s. of ( |3.34j ) 
by 3> e (x), we have to prove the relations: 

(a) lim £ _ T0 $ £ (x) = $±(x); 

(b) lim £ ^ T0 E xe z d IM*) - $±(x)| 2 = 0. 



where $±(x) are defined in ( 3.16 ). Both facts are proved in the Lemma [TT] below. Theorem 
373] is proved. 

Remarks. 1). Functions ^±(x, k) are analogs of the Sommerfeld solutions, which appear in 
the scattering theory for potentials decaying in all directions and which provide a complete set 
of generalized eigenfunctions for the part of the spectrum that coincides with the spectrum of 
the Laplacian [^, |2^|. Likewise, ( p. 15 ) is an analogue of the Lippmann-Schwinger equation of 



scattering theory. 

2). According to formula (|3.13|) , \l/±(x, k) depends on the component x 2 of x = (xi,x 2 ), X\ G 
Z dl , x 2 G Z d2 via the product of e lk2 ' X2 and of a 1-periodic function of the argument a-x 2 , i.e. of 
a quasi periodic function of x 2 G r L d2 (recall that we assume in this section that the components 



of the vector a in (|1.4j ) are rationally independent). This fact is in agreement with the widely 
accepted idea, according to which generalized eigenfunctions of absolutely continuous spectrum 
of differential and finite difference operators with almost periodic coefficients have the "almost 
Bloch" form, i.e the form of the product of a plane wave and an almost periodic function with 
the same frequencies as the coefficients (see e.g [p2|). 



3). According to formula (|2.18|) , if \E\ > u, the Green function Gq ' (x; £ + iO) of the v- 
dimensional Laplacian decays exponentially and if \E\ < v it decays as l/j^l^ - for v > 2 
(in the one dimensional case for \E\ < 1, G^\x; E + iO) behaves as e^^*"^, where r](E) is 
a real valued function, see formula (|3.37| ) and ( |3.38| ) below ). As m varies the expression 



Ed(k) — Ed 2 {k 2 — ma) has values inside (— di, d\) as well as outside this interval, then the Green 
function 

entering the expression (|3.13|) , may be exponentially decaying or slowly decaying (i.e. as 
l/lxj 1 ^ - ). In other words we can write, say for 

*(x, k) = e 2 ^ k - x + ^(x, k) + ^ sur/ (x, k), (3.35) 
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where ^/ vo i is the part of the sum in ( |3.13| ), containing only slow decaying terms, and ^ sur f is 
the part, containing the exponentially decaying terms. 

Recall now the definition of the surface states according to [[HJ (for other definitions see 
1,0,0)- 



Definition 3.1. Let ^> e be a generalized eigenf unction ^e, corresponding to a point E of the 
spectrum of the operator H of (\1.1\) - t \1.3j) . We say that \1/ E is a surface state, if for any e > 
we have 

sup (l + lzal*/ 3 ^)- 1 \Ve{(xi,x 2 ))\ 2 < oo. (3.36) 

x » eZda z ie Z<«i 

Since the part e 2j7rk x + ^/^(x, k) of the solution ( |3.35| ) ^(x, k) is not decaying in the X\- 
direction, the solution is not a "surface" state but a "volume" state. Hence, we can say that 
Theorem |3]3| above implies the existence of the volume states for all E G (— d, d). Theorem |3.4| 
below implies that these generalized eigenfunctions are complete in the interval (—d,d). We 
conclude that there is no surface states in the interval (—d, d) of the spectrum of the operator 
H in the considered case of quasi periodic surface potential (|1.3|) - (|1.4j ). However, despite 
that surface states are absent, the volume states (|3.35|) contain both a term, e 2l7rk ' x + ^ vol ('x, k) 
which slowly decays or even only oscillates in \xi\ , and a term, 1 i r sur /(x, k), which exponentially 
decays in |xi|. They are respectively the superposition of reflected or transmitted waves which 
propagate inside the bulk and of waves which propagate only along the subspace Z d2 . 

4). The scattering interpretation ( |3.35| ) of generalized eigenfunction ( |3.13| ) allows us to intro- 
duce transmission and reflection amplitudes and coefficients (the latter as square of the modulus 
of the former). Consider the simplest case of d\ = 1 and recall that: 

• p iri(z)\x\ P iri{z)\x\ 

Gq \x; z) = ——- = - -=== , (3.37) 

Sill T](z) y/z 2 -l 

where — cos rj = z, or 



Tj(z) = -i \og(-z + Vz 2 - 1), (3.38) 

and we use the branch of the logarithm that has the cut along the negative semi-axis and the 

branch of \J z 2 — 1 fixed by the condition \/ z 2 — 1 = z(l + 0(z~ l )),z — > oo. In particular 
^sr)(z) > for > and 

(0,vr), \E\<1, 

r]{E + iO)e{ Tr + iR+,E>l, (3.39) 
+iR+, E < -1. 



Combining these formulas and ( |3.13|) , we can present ^ vo i (x, k) in (|3.35| ) for d\ = 1 as 



* w i(x, k) = ^ m (k)e ir?m(k)|a;i|+2m(fc2 ^ m) - X2 , (3.40) 

m 

where denotes the sum of those terms in ( |3.13|) for which r] m (k) := r/(A m (k) + iO) is real, 
and A m (k) is defined by the equation: A m (k) = Ed(k) — ^-1(^2 — ma). Recall that in the 
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one-dimensional scattering problem for the potential vS(x),x G Z, the Sommerfeld solutions 
are (cf (f£50p , and (ggg )): 



if + sin27rfc 

where 77 (A;) = ^(cos27r/c + zO), G T. Hence in this case 

... sin27rA; .,. if 

t(fc) = . - r _(k) 



iv + sin 2nk iv + sin 27r/c 

are the transmission and the reflection amplitudes. This makes natural to view 

t = l + #o(k), r = $ (k), 

where \l/o(k) is given by (|3.40| ), as the transmission and the reflection amplitudes of the plane 
waves scattered by the surface potential ( |1.3j ) and propagating in direction k of the incident 
wave and in the opposite direction. Likewise it is natural to view the coefficients \I/ m (k),m > 
1 of (|3.40|) as the transmission and the reflection amplitudes of the scattered plane waves 
propagating in the directions (r] m , k 2 + ma) and (—r] m , k 2 + ma) respectively to the right and 
to the left of the plane x± = 0. This scattering theory interpretation of the solutions (|3.13| ) is 



in agreement with the form of the scattering matrix S in our case. We use the general formula 



(see [H, formula (4.2.30)): 

-2m) s - lim / 5 £2 {H - A)T(A + ie 1 )S Q {dE), (3.41) 



s = 1 - r, r 

£2->0 ' 



where S e (A) = (2m)~ 1 [(A + ie)' 1 - (A- is)' 1 }, T(z) is defined in (|J), £ is the resolution of 
identity of H , and the limits have to be carried out in the following order: first E\ — > 0, second 
£2 - > 0. Formula ( |3.41| ) implies that for any sufficiently smooth function / on T d we have: 

(T/)(k) = (-2m) s- lim / 5 £2 (E d (k) - E d (p))T(k,p;E d (p) +ie 1 )f(p)dp. 

By using formula ( |2.29| ) for the kernel of the T-operator, it can be shown that the generalized 
kernel T(k, p) of the T-matrix of Q3.41| ) is: 

T(k,p) = -2m8{E d {k)-E d {p))T{k,p;E d {p)+iO) (3.42) 

00 

= -2m5{E d (k) - E d {p)) t m {h; E d {k) + i0)5(k 2 + ma - p 2 ). 

m=0 

Now we formulate and prove the lemma that was used in the proofs of assertions (iii) and (iv) 



of Theorem ET3 



Lemma 3.1. Let $x,2 : T dl ' 2 — > C be smooth functions. Assume that the support of $1 does 
not contain the critical points of E dl : 

supp $! f|{k! G T dl : VxE dl (kx) = 0} = 0. (3.43) 
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Set for e ^ 0: 



*£ d (k)+ l£ (x,k)<l(k)rfk, 



where ^ z (x, k) is defined in $TQ) , t \3.10j ), and in Q3.19{ ), and $(k) is of the form $((fci, k 2 )) 
$(fci)$(fc 2 ). T/ien: 



|$ e (x) -$±(x)| 



0. 



(3.44) 



where $±(x, k) are defined in ( $.12j and in ( \3.16] ). 
Proof. By using ( |3.19| ), we find that for any e 7^ 0: 



$,(x1 



$(x) 



00 

E 

m=0 



T d l 



dpi$i(pi) / dp 2 e 2i7r(fc2 - ma):C2 <i) 2 (p 2 



T d 2 



(3.45) 



x t m (/c 2 ,£ d (k) + ie)G[ ) dl) (£ 1 ,£ d (k) - E d2 (k 2 -ma) + ie) 



where $ is the Fourier transform of $. The integrals and the series in this formula converge 
and can be written in any order because of the bound (|3.3|) applicable in view of ( |3.43|) . The 



integral representation (|2.18| ) for G Q dl ^ allows us to rewrite the last formula as follows: 



$(x) + 



„2i7rk-x 



^ £ (k)rfk, 



where 



*«(k) 



E 



dpi 



$l(Pl)$2(&2 +H 



— j Tdl £ d (k) - (pi) - £ da (fc 2 + ma) + 
x t m (k 2 ,E dl (p 1 ) + E d2 (k 2 + ma) +ie). 



(3.46) 



(3.47) 



This series and the integral are convergent because the modulus of the denominator is bounded 
from below for e 7^ 0, and because of bound ( p.3|) . 

Now we will prove that for any k G T d , the limits lim^zpo ^e(k) = ^(k) exist and that the 
convergence is bounded. Consider the case for the sake of definiteness. The building block of 
the coefficient t m (k 2 , E + ie) in ( ET37D is the function 70(^2,^ + ^) = G[, dl) (0, £?- E d , 2 (k 2 ) +ie)). 
This function is real analytic in k 2 G TP 2 (see (|3.11|) for the definition of TP 2 ), and in E G 
(— d + 7, d — 7) for any fixed (small) 7 > 0, (see ( |2.18| ) and ( |2.19| )). By using identity ( |5.8j ) for 
(?Q dl \ we can write the mth term of formula ( |3.47| ) as: 



ma 



POO 

/ dte~ £t ~ it( - Ed( - ]i) ~ Ed 2 < - k2+ma ^& 2 (k 2 + 
Jo 

dpA l { Pl )e- itEd ^ pl k m {k 2 ,E dl {p 1 )+E d2 {k 2 + ma)+te)). 



(3.48) 



Since the support of $1 does not contain critical points of E dl and since G m (k 2 , E + ie) is real 
analytic in k 2 G T^ 2 and in E G {—d + 7, d — 7), 7 > for all £ > 0, we can integrate by parts 
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twice in components of p\ G T dl , and obtain an expression of the form t _2 $ m (k,p 1 , e), where 
$ m is bounded in k G T d , p x G T dl and e > 0. This allows us to make the limit e — > +0 in 
(|3.48f) and obtain a bounded in k expression. 



Besides, $ m is a linear combination of the first and second partial derivatives in components 
of pi G T dl of the integrand in ( 3.48|) . The derivatives are linear combination of products 
of bounded (and smooth) in k G T ,p% G T dl for e > 0, and independent of m func- 
tions, multiplied by the first and the second partial derivatives in components of p\ G T dl of 
ni^i b(k 2 + la, E dl {p 1 ) + E d2 (k2+ma)+ie)). This leads to the bound |$ m (k,p 1 ,e)| < cim 2 e~ C2m 
where Ci, c 2 > and are independent of m, k,pi and e. The bound allows us to make the limit 
e -> +0 in (0T47|) for any k G T d : 



lim * e (k) 

e— >+0 



¥_(k), 



and to obtain the bound |\l/ e (k)| < const, valid for any e > and k G T d . Now the Lebesgue 
dominated convergence theorem and relation (|3.12|) proved above lead to the representation: 



#_(x) = $(x) + / < 



2i7rk.x 



tf_(k)dk. 



(3.49) 



Subtracting this relation from ( |3.46 ) and applying to the result the Parseval equality, we obtain 
that: 



*-(x)l 2 



Jd 



|i e (k) -^_(k)| 2 dk. 



(3.50) 



Thus ( |3.49| and the Lebesgue theorem imply ( |3.44j ). Lemma is proved. 



Theorem 3.4. Let H = Hq + V be the self-adjoint operator on l 2 {7j d ), defined by ( \1.J\) - 
( \1.3j ) in which the vector a G M^ 2 has rationally independent components. Then the family 
{\E^(x,k);x G Z d } k( -f d , defined in Theorem (see flSlfy, ( \3.12j ), and ( \3.13j) ), is the complete 
system of generalized eigenf unctions of H in the part (—d,d) of the spectrum of H, i.e.: 



(i) for any f G I 2 {If), the series 



F±(k) = £* ± (x,k)/(x) 



(3.51) 



converges in Z 2 (Z^ 



(it) if £h{A) is the spectral projection of H , corresponding to the closed interval A = [a,b] C 
(— d, d), then 



l^(A)/f= / \F ± (k)\ 2 dk; 
Jket d :E d (U)eA 



(3.52) 



where Ed(k) is defined in ( \2.1Q )); 
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(Hi) the following relation is valid 

\\H£ H {A)f\\ 2 = [ \E d (k)F ± (k)\ 2 dk. (3.53) 

JkeT d :E d (lc)eA 

Proof. We write the Hilbert identity for the Green function G(x, y; -21,2), 3^1,2 7^ 0: 

G(x, y, zi) - G(x, y, z 2 ) = {z x - z 2 ) G(x, s; «i)G(y, s; z 2 )- (3.54) 



By using the Parseval equality for the Fourier transform with respect to the variable s in the 
r.h.s. of this identity, we rewrite it as follows: 



/ dkG(x, k, zi)(7(y,k, z 2 ) 



where G(x, k, z) is the Fourier transform of G(x, y, z) in the second variable y, defined in ( |3.9|) . 
Multiply now resulting relation by /(x)/(y), where / has compact support in Z d and sum over 
x, y G Z d . This yields: 

»' - G( - ,)/ ' /) = l V)-^ W -^ , 

where 

F*(k) = £tt 2 (x,k)/(x), (3.55) 

xez d 

and \l/ 2 (x, k) is defined in (|3.10| ). Setting zi = z 2 = E + ie, e > 0, we get: 
1 



Q(G(E + ie)f,f) = - dk 2 \F E+ UV\ 2 . (3.56) 

7T 7T Jfd (-^d(k) - E) 1 + £ 2 



A = [a, 6] G (— d, d), we obtain in the l.h.s. of the resulting relation the expression ||£f/(A)/|| 2 . 
can be continued in z to the real z = Ed(k) + iO G A, and that the continued function is 
uniformly continuous in k G {k G T d : Ed(k) G A}, where T d is defined in ( ft.ll ). Since / 



is of compact support in Z d , it suffices to show that ^(x, k) possess this property for any 
fixed x G 7L d . But this fact is proved Theorem |3.3| . Thus we have established (ii) for the case 
where /'s of finite support. The extension to /'s belonging to / 2 (Z d ) is based on the standard 
arguments of spectral theory (see e.g. |25|, |23|). This proves assertions (i) and (ii). As for 



assertion (iii), it follows from (ii) and from the spectral theorem. 



4 The Periodic Case 

In this section we consider the operator H = Hq + V of fll.l]) - (|1 . 5|) in which d\ = d 2 = 1 and 
a is a rational number: a = p/q, p G Z, q G Z \ {0}, i.e. for periodic potentials v of fll.4p . 
We show that in this case the whole spectrum of H is absolutely continuous and we construct 
corresponding generalized eigenfunctions. It turns out that there are two types of generalized 



22 



eigenfunctions. Both types have the Bloch-Floquet form in the longitudinal coordinates x 2 but 
behave differently in the transverse coordinate x\. 

We will follow the same strategy as in the preceding section namely the construction of gener- 
alized eigenfunctions based on the formulas for the Green function of Section 2 and on formulas 
(|3.9| ), ( 3.1C ) and ( |3.12| ) of Theorem |3.3| . Thus we have to analyze the behavior of the Green 
function as the spectral parameter tends to the real axis. Our first goal is to find the set of 
energies for which the limit G(x, y, E + iO) exists and is bounded, i.e. the purely absolutely 
continuous part of the spectrum. We shall see that unlike the quasiperiodic case, where this set 
is [—d, d] , in the periodic case the whole spectrum is pure absolutely continuous. The spectrum 
which lies outside [—d, d] consists of surface states only. As for the part in the interior of [—d, d] , 
it consists of the volume states whose energies occupy the whole interval [— d, d], and of the 
surface states that may exist under certain conditions. 

For any z6C, Qz 7^ 0, and m — 1, q define the function: 



P m (k 2 ; z) = a m Yl b(k 2 + la; z),Wk 2 e T, 



(4.1) 



1=1 



where a and b are defined by ( P-llQ , (|2.12|) , (|2.25 ) and T = (0, 1]. Then, by using Lemma [T2 
we obtain for a = p/q: 



((foovxAfe) = <? q n ^ + ia ~> z )^) = p ^ z Mk 2 ), 



1=1 



(4.2) 



where the operator u is defined in ([2.9Q. We conclude that (bu) q is a multiplication operator 
by the function P q in the space L 2 (T). 

Theorem 4.1. Let H = H + V be the operator defined by - ( \l-o] ) in which d\ = d 2 = \ 

and a=p/q,p^Z,q^'Ij \ {0} is a rational parameter. Then the Green function G(x, y; z) = 
(H — 2) _1 (x, y), x, y G 1? of H can be written in the form: 



G(x,y;z) 



Gf(x 



Gq(xx', z + cos27r/c2)G , Q i ' ) (yi, z + cos27r(/c 2 + ma))e 



y;z) + V / dk 2 e 2 ™ k ^~y*H m (k 2] z) 

m=0 



(4.3) 



-2iirmay2 



where 



tm(k 2 ; z) 



9 



g%(k 2 ;z) + i 
( -1, 

1 



(4.4) 



2ia 



1 - P q (k 2 ; z) g^foih + a; z) + 1 
1 2ia 

k 1 - P q (k 2 ; z) g^o(k 2 + ma; z) + i 



m = 0; 
m = 1; 



P m -i(k 2 ;z), m > 2, 



■ \z), V 



1, 2 is the Green function (\3.31Q of the v- dimensional discrete Laplacian, %(.] 



and b(.; z) are defined in ( \2.23j) and in (\2.25 ). 
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The proof of the theorem is based on the same argument as that used in the proof of Theorem 

o 

Formulas (|4.3|) and (|4.4j) suggest that the spectrum a(H) of H contains the set S = {E G R : 
3k 2 G T;P q (k 2 ,E) = 1}. We prove below that indeed, the limit G(x, y, E + iO) exists and is 
bounded for all E G S \ D where D is a discrete set. 

For any 1>7>0, G R and I = 1, q, define the sets: 

Kl f (E)={k 2 eT: E + cos2ir(k 2 + la) G [-1 + 7 ,1 -7]}, (4.5) 

and 

K y (E) = \jK! r (E), K;(E)=T\K,(E). (4.6) 
1=1 

It follows from formula (|4.2j) , Lemma |5l|, and from the argument of the proof of Theorem [3.1| . 
that for any < 7 < 1 there exists #(7) > such that the inequality sup e>0 |P g (fc 2 , E+ie)\ < 1—5 
is valid uniformly in k 2 G K 1 (E). This means that the function (1 — P q (k 2 ,E + iO))' 1 is well 
defined and bounded on the sets, 

K{E)= |J K y (E) 

0<7<1 

and that possible singularities of this function which are given by the "band-equation": 

P q {k 3 ,E) = l, (4.7) 

where P q (k 2 , E) = P q (k 2 , E + iO), are localized on K C (E). It is natural to think that energies, 
satisfying the band equation (|4.7| ) for some k 2 G T belong to the spectrum of H . The following 
proposition describes properties of solutions of the band equation. 

Proposition 4.1. For any 2 < q < 00 the band equation (\4-7\) admits a finite number N' of 
positive solutions < E 1 (k 2 ) < ... < E N i q {k 2 ) < 00 (the positive energy band functions), and 
a finite number N q of negative solutions —00 < E_ N ^i(k 2 ) < ... < E_i{k 2 ) < (the negative 
energy band functions) . 

The functions Ej, j = —N q , ... — 1, 1, ...A^ are 1/q-periodic in k 2 , and are real analytic in the 
interior of their respective domains T>j C T (each domain T>j is a closed subset off). 

Moreover, the band functions are separated in the sense that: 

(i) for any j = —N", ... — 1, 1, ...N q there exists a finite subset Vj of T>j, such that for all 
k G T x (Vj\V-) we have: 

\E j (k 2 )-E 2 (k)\>0, (4.8) 
where E 2 (k) = — cos27rA;i — cos27rA;2/ 

(ii) there exists a positive constant r\ q > such that for any j,j' = —N q , ... — 1, 1, ...N' q , j 7^ j' 
we have: 

inf \E j {k 2 )-E jl {k 2 )\>r ]q >Q. (4.9) 
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The proof of the proposition will be given after the proof of Theorem WM. 



The band function Ej, j = —N{j,... — 1,1,. ..N^ defined in Proposition ( f4.1| ) determine the 
band-gap structure of the spectrum of the periodic in x 2 operator H in the following sense, 



Theorem 4.2. Let H = Hq + V be the operator defined in Theorem Then for all rational 
parameter a, g ^ 0, and us G [0, 1] the spectrum c(H) of H is a finite union of closed intervals 
(energy bands): 



N' 

<i 



a(H)= [j Ran Ej U [-2,2]. (4.10) 



-Nil 



The assertion that (—d, d) is in the spectrum of if is a consequence of Theorem |2.2|, the rest of 



the theorem will be proved after the proof of Theorem [4J) 
Let us now define the set E c of critical energies as 

dE 

S c = {E G R : 3j G {-N" N'}, 3k 2 e T, Ej{k 2 ) = E and either -77^(^2) = or k 2 e dVj}. 

dk 2 

Denote 

D = £ c \J{-d,d}, (4.11) 



and notice that because of Proposition |4.1| D is a discrete subset of R. 

Theorem 4.3. Let H = Hq + V be the operator defined in Theorem and let G(x, y; z) be 
its Green function. Then for any rational a, g G R, and u G [0, 1) the limit G(x, y; E + iO) 
exists and is bounded for any E G cr(H) \ D and x, y G 1? , where D is defined in (\4-ll ). In 
particular the spectrum of H is absolutely continuous. 

Proof. For any E G o~(H) \ D set z = E + ie, e > and fix < 7 < 1. By using formula (|4.3|) 
we can write that 

G'(x,y;z) = G' 1)7 (x,y;z) + G 2i7 (x,y;z), x,yGZ 2 , (4.12) 

where 

r q 

G li7 (x, y; z) = 4 2) (x -y;z)+ dk 2 e 2 ™ k ^~y^ V t m {k 2 - z) (4.13) 

J K r/ (E) m=Q 

x G$\xi; z + cos27rA; 2 )G'J 1) (yi; z + cos2tt(A; 2 + ma))e" 2! ™ s , 

Kj(E) and t m are defined in (|4.6| ) and in ( |4.4|) , and 

G 2 , 7 (x, y; z) = G(x, y; z) - G?i, 7 (x, y; z). (4.14) 

Since the inequality sup e>0 \P q {k 2l E + ie)\ < 1 — 5 is valid uniformly on K y (E), the same 
arguments as in the proof of Theorem |373] imply that the limit Gi i7 (x, y; E + iO) exists and is 
bounded. 
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Hence, to prove the theorem we have to show the same property for the term Cr2, 7 (x, y; z) of 
(|4.12 ). We first note that by Proposition [4 . 1| this term can be rewritten as 

n i \ f ji #2,7 (x, y, k 2 ; z) . . 

G 2 , 7 (x,y;z) = / dk 2 ,y -, (4.15) 

where for any < 7 < 1, e>0 and (x, y) G Z 2 x Z 2 , <72, 7 (x, y, ; z) are smooth functions on 
K^(E). Now in order to compute the integral in the r.h.s. of fl4.15|) , consider the level sets: 

S 3 = ^(£,7) = {k 2 E K°{E) : E 3 (k 2 ) = E}, j = -N' q ',... ,N' q , 

and the following neighborhoods Uj of Sj: 

Vj = = {h e K C JE) : \Ej{k 2 ) - E\ < V }. 



If 77 is small enough, then Proposition |4J] implies the relation: Uj H vy = if j 7^ j'. Thus to 
prove that G 2 ^(x, y; E + is) exists and is bounded as e — > 0, it suffices to show that this holds 
for 

n ( tp , • \ f ai #2, 7 (x,y, k 2 ;E + ie) . 

G 2iltj (x,y;E + te) = / dfc 2 — ^ ^-7 — r , 3 = -N q ,...,N q . (4.16) 

Since 77 is small enough and E D,we can parameterize z/j by the local coordinate E defined by 
the relation E = Ej(k 2 ). Denoting ipj the respective change of variables and J Vj its Jacobian, 
we have 

G 2 , 7J (x, y-E + is) = [ V dE ) j . = ^ . . . j # (41?) 

J-V l - P q°VA E ) 

Suppose now that 77 and e are so small that we can write: 

1 - P q (<pj {E),E + ie) = {E-E-ie) Pj {E;E + ie),E e [-77,77], 

where pj, j = —N' 1 ,--- ,N' are smooth and non vanishing functions on the interval [—77,77] 
such that 

\ Pj (.,E)\>C\d E P q (.,E)\+0( V ) + 0(e) 

for some strictly positive constant C . Moreover it follows from the proof of Proposition [11] (see 
formula |4.45| ) that 

\d E P q {^{E),E)\^^ £€[-77, 77]. 

Then standard arguments imply the existence and the boundedness of Gr2, 7 j(x, y; E + iO), 
j = —N q , ■ ■ ■ ,N q hence the existence and the boundedness of Cr2, 7 (x, y; E + iO). The theorem 
is proved. 

The last theorem together with the arguments of the proof of Theorem |3.2| lead to: 
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Theorem 4.4. Under the conditions of the Theorem [O, the wave operators Sl± for the pair 
(H, Ho) defined in - with a rational a exist and are complete for any closed interval 



[a, 6] c (-d, d) \ U 



N' 

■jl- N » Ran^- 



Our next theorem shows that surface states (see definition ^.1| ) exist and are bounded. They 
can be labelled by the " quasi- momentum" k 2 6 T/g, such that respective eigenvalues are given 
by the band functions: E = Ej(k 2 ). The "volume" states that do not belong to Z 2 (Z) in x\ are 
labelled by the "momentum" k e T 2 , such that the corresponding eigenvalues are given by the 
dispersion law of the Laplacian: E = ^(k) = — (cos27r/ci + cos2-7rA;2)- We consider here only 
the non-degenerate case, i.e. the case where chosen pairs (k 2 , E = Ej(k 2 )), and (k, E = E 2 (k)) 



are such that Ej(k 2 ) 7^ -^(k). By Proposition ^j] this property is valid for all energies except 
a finite set. 

Consider the set: 



T 



{k=(k 1 ,k 2 )et 2 , fejGT.], j 



-N" 

1 v q 1 •••) 



N' 



where T 2 is defined in (|3.11| ), Tj = T>j \ Vj, and T>j,T>'j are defined in Proposition (|4.1|) , and 
the set 



N', 



T 



|J {k = (k 1 ,k 2 )ei 2 ,k 2 eJ\V' :j }. 



Hence the set of degenerate energies is 

£' c = {Ee E, 3k = {k u k 2 ) e T 2 , 3j 



-NZ,...,NL,E = E 2 OL) = E j (k 2 )}, 



By Proposition fl4.1|) £' c is a discrete set as well as the set 

D' = DU E' 



where D is defined in ( [4.11 



Theorem 4.5. Let H = Ho + V be the operator defined in Theorem G(x, y; z) be its Green 
function, and G(x, k; z) be defined in ( \3. Q ). Then: 



(i) for z = ^(k) =F is the limits 



*«,±(x,k) = lim #*(x,k) 



= lim ±zeG(x, k; (E 2 (k) =f ie)), (4.18) 



exisi /or a// k 6 T 2 , are bounded in x e Z 2 /or any k 6 T 2 , are continuous in k on any 
compact subset 0/T 2 /or any x G Z 2 ; and satisfy the Schrddinger equation: 



({H + V0*„, ± )(x, k) = £ 2 (k)tt„ i± (x, k); 



(4.19) 
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(ii) For z = Ej(fa) =F is, j = —N 1 ' ...,N' the limits 



*«j,±(x, fa) = lim Tielj(k; z)G(x, k; (^(fc 2 ) =F ie)), 



in which 



n l/2 

/,•(]*; z) = (£ 2 (k)-z) 



t A |(i&(k)-s)p» 

exisi /or any k = (fa, fa) G T 2 ; are bounded in x G Z 2 /or any fc 2 G Tj ; are continuous in 
fa on any compact subset of Tj and satisfy the Schrddinger equation: 

((H + V)V s , jj± )(x, fa) = %(A; 2 )* Sli ,±(x, fc 2 ). (4.20) 
^'mJ ^ r s ,i,±( - ) fe)) ^2 ^ are surface states in the sense of Definition 3.1. 

N' 

Remarks. 1). It can be shown that for all k G T 2 such that -E 2 (k) G (—d,d) \ [Jjl_ N „ R&nEj 
the function *& Vj ±, defined by (|4.18|) , is the unique solution of the integral equation: 

*(x, k) = e 2 ^ k - x - G ? ] (x - y; E 2 (k) t tO)V(y)^(y, k), (4.21) 



that has to be understood in the same way as in Theorem |3.3| fii). On the other hand, it is easy 
to check that for any j = —N", N' q , k = (fa, fa) G T 2 and Ej(fa) G" [—d, d], \1/ s j i ±(x, fa) is a 
solution of the homogeneous integral equation: 

tt(x, fa) = - G o\* ~ y; Ej(fa) t iO)V(y)*(y, fa). (4.22) 

2). One can view the above results from the point of view of the direct integral decomposi- 
tion technique for finite difference operators with periodic coefficients Namely by using the 
periodicity in x 2 of the operator H with a = p/q, we can write the direct integral decomposition 



H= / H q (fa)dfa. (4.23) 

J\k 2 -l/2\<l/2q 

Here H q (fa) is the selfadjoint operator defined by the restriction of H to the linear manifold 
of functions ^/fc 2 ( x ) = e 2i7rfc2I2 $fc 2 (x), where $fc 2 is g-periodic in x 2 . Thus H q (k 2 ) acts in the 
strip {%i G Z, 1 < x 2 < q}, and is the perturbation of the respective Laplacian by the q 
rank potential (|1.3|) with 1 < x 2 < q. This implies that the spectrum of H q (k 2 ) consists 
of two parts. The first is the absolutely continuous component: the union of values of the 
functions -cos2vr£;i - cos2vr(A; 2 + l/q), fa G T, 1 = 1, ...,q and fa G [1/2 - l/2g, 1/2 + l/2g) 
is fixed, the corresponding eigenf unctions are deformed plane waves in x±. The second part is 
discrete spectrum, consisting of N q < q eigenvalues Ej(fa), lying outside of the above absolutely 
continuous spectrum, and having exponentially decaying in x± eigenfunctions.As fa varies in 
the direct integral the absolutely continuous spectrum of H q (k 2 ) gives rise to the volume states 
of the operator H, while the discrete spectrum of H q (k 2 ) gives rise to the surface states. 
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Proof of Theorem g73f Take (22, k) G a(H) x f 2 , z = E ± ie,e > and denote *«(x,k) = 
(22 — z)G?(x, k; 2), where G(x, k; z) is defined in (|3.9|) . We know from the proof of Theorem |3.3| 
that if for any x6Z 2 the limit \I/e(x, k) = lim £ ^o ^(x, k) exists, then ~§>e is a solution of the 
Schrodinger equation 2f\l/£ = E^g. By Theorem |4.1| we can write the representation: 



G(x, k, z) 



E 2 (k)-zl 



2j7rk-x 



(4.24) 



+ t m {k 2 - ma; z)G ( q ) (an; z + cos 2tt(A;2 - mtt))e 2ff(trma)l2 



m=0 



Choose first a pair (k e T 2 , £ = ^(k)), as it was done in the proof of Theorem |3]^ for the 
quasiperiodic case. By Proposition [4.1| the denominator 1 — P q (E, k 2 ) in t m of ( |4.4| ) is nonzero 
and we obtain from ( |4.24| ): 



*«,±(x, k) 



2i7rk.x 



+ ^m(^2 — wia; z)Gq \xi\ z + cos27r(A; 2 — ma)) 



m=0 



(4.25) 

2in(k2— ma)X2 



2=B 2 (k)T*0 



This proves the first assertion of the theorem. 

Consider now the case where k = (ki, k 2 ) G T|, and 22 = Ej(k 2 ) for some j = — iV^ , N' . We 
know that the pair (k, E) is such that E + cos27r(/c2 + lot) G" (—1, 1) for any / = 1, q. Hence, 
by using the separability property (|4.9|) and the periodicity of the 22j's, given by Proposition 
P~T| , we find that E + cos27r(&;2 + la) G" [—1, 1] for any I = 1, ...,q, i.e. all that these energies 
belong to the resolvent set of the 1-dimensional Laplacian. 

This observation implies the existence of the limit 



* a j,±(x, k 2 ) = hm Q [J(k, z)V K (x, k 



provided that the limit 



lim 



e^+o I - P q (k 2 , Ej(k 2 ) t is)' 
exists. This can be proved by using the relations 

1 - P q (k 2 , E 3 {k 2 ) t is) = ±ied E P q {k 2 , Ej(k 2 )) + 0(e 2 ), 

valid for sufficiently small e, and the relation d^Pq 7^ 0. Now, it easy to verify that 

# sJ ,±(x, k 2 ) = Ijih) tm(k 2 - ma)G ( i\x l] E,{k 2 ) + cos2vr(A; 2 - m«))e 2w(fe -™ )l2 , 

where 



(4.26) 



(4.27) 



m=l 



dk\ 



1 



\(E 2 (k)-z) 



1/2 



(4.28) 



Z=Ej(k 2 )^iO 
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. -1 

ti(fc 2 ) = 2icx( (g%(k 2 ; z) + i)d E P g (k 2 , z)(g%(k 2 + a; z) + %)) 



z=E j (k 2 )^iO 



and for m > 2 

t m (k 2 ;E) = 2iga(ig%(k 2 ;z)+i)d E P q (k 2 ,z)(g%(k 2 + ma;z)+i)\ (4.29) 
x P m _i(k 2 ;z) 



z=E ] {k 2 )TiO 



By using the same argument as that in the proof of ( 4.19|) , we find that ^ s ,j,± satisfies ( |4.20| ) . 
Let us prove now that *& s ,j,± is a surface state. We know that E + cos2i[{k 2 + la) G" [—1, 1] for 
any I — 1, • • • , q. Since all these energies are in the resolvent set of the 1-dimensional Laplacian, 
each term of the sum of the r.h.s. of ( }4.27|) decays exponentially with respect the transverse 
coordinate x% G Z. Since the number of these terms is finite, we conclude that for any x 2 G Z, 
^,±(■^2) G Z 2 (Z). The proof of the theorem is complete. 

We can now use the last theorem, where we have constructed the generalized eigenfunctions 
QOp and to prove Theorem [O. 



Proof of Theorem [4.2| . It follows from the proof of the Theorem |4.3| that er(if) c [— d, d] U 



TV' 

(Uj=-jv" R- an Ej)- Hence we have to prove the opposite inclusion. For the part [—d,d] of the 

spectrum the inclusion was proved in Theorem |2T2| . So assume that E G U^-iv, Rani^j \ [— d, d] 
is such that there exists a surface state \P s (x) satisfying the Schrodinger equation: (ff\I/ s )(x) = 
_E\I/ S (x). We apply again the H. Weyl criterion, setting 

* n (x) = l n (:r 2 )^ s (x)/iV n ; N n = ||l n #J, 3(za) , 

where l r is the indicator of the ball {x 2 G Z : \x 2 \ < n}. A straightforward calculation shows 
that C\n 1 / 2 < N n < C 2 n 1 ^ 2 as n — > 00 for some strictly positive constants Ci j2 , and that 

E^ n (x), \x 2 \ < n; 
(HV n )(x) = { A n {y), n<\x 2 \<n + l; 

£tf„(x) = 0, |x 2 | > n + 2, 

where ||^4^||z 2 (s 2 ) = Ofn^ 1 / 2 ), n — > 00. It is easy to check that \l/ n is a Weyl sequence for i7 
at the energy i?. This proves the theorem. 

Our next result concerns the completeness of the system of generalized eigenvectors ([4.25 ) and 
( 4.27p ) defined in Theorem | 



Theorem 4.6. Let H = H Q + V be the selfadjoint operator in / 2 (Z 2 ) defined in Theorem 

N' 

Consider the family T = {^(x, k), x G Z 2 } kg ^ 2 Uj=-jv"{^»j( x j ^ 2 )> x e ^l^eTj? defined by 
( U-^3j ) and by ( \4-27{) . Then J 7 is a complete system of generalized eigenfunctions of H in any 
sufficiently small interval A of a(H) such that A n D' = 0, i.e. 

( i) for any f G Z 2 (Z 2 ) the series 

F„(k)= ^^(x,k)/(x), kef 2 , 
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and 

F sJ (k 2 ) = ^ SJ -(x,fc 2 )/(x), h G Tj, j = -iV;, N> 

x£Z d 

converge in Z 2 (Z 2 )/ 

(ii) if 5ff(A) zs £/ie spectral projection of H corresponding to the interval A G o~(H), then 

N'„ 

\£ H (A)f\f 



[ |F,(k)| 2 rfk+ V / \F hS (k 2 )\ 2 dk 2 ; 

J{kef2 :j B 2 (k)eA} n=-N" J{k2&Tj:Ej(k2)eA} 

J 1 V 



fmj /or t/ie same interval we have 

\\H£ H (A)f\\ 2 = [ \E 2 (k)\ 2 \F v (k)\ 2 dk 
J{kef 2 :E 2 (k)eA} 

N' 



+ J] I \E,(k 2 )\ 2 \F SJ (k 2 )\ 2 dk 2 . 



Proof. For any compact interval A C cr(if) \ £)' consider the sets: 

1/ = {k G f 2 : E 2 (k) G A}, ^ = {k = (ki, k 2 ) G T 2 : Ej(k 2 ) G A}, j = -JVj', N' q . (4.30) 

Proposition f4.1| implies that there exists a constant rj > such that: 

min inf |£ 2 (k) - EAk 2 )\, min inf |E 2 (k) - EAk 2 )\ > rj. 
j kgi/nT 2 j ke^nf 2 

Notice that ?7 depends only on the dist(A, D'). Moreover, if A is sufficiently small, then the 
sets = —N",...,N' are disjoint. The subsequent argument uses this property of v, and 

We will follow now the proof of Theorem |3.4| , Hence we have to prove assertion (ii) first for a 
function / with compact support. We have for z = E + ie, where E G A and e > 0: 

isKcw/,/) = ijfdk |(£ ,^_ t)| , |F.(k)|' 

where 

F,(k) = ]T *,(x,k)/(x), *,(x, k) = (£ - z)G(x, k; z), (4.32) 



E = E 2 (k) or E = Ej(k 2 ) and G(x, k; 2) is defined by (gj) and (Ml)-(M)- Since for every 
k = (ki,k 2 ) G ^ and £ 6 A, P q (k 2 ,E) — 1 is not zero, the limit lim e ^o ^s+ie(x, k) exists for 
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any x£ Z 2 uniformly in k = (fci, fca) £ ^ and in E G A. Applying to the first term of the r.h.s 
of ( |4.31|) the operation lim £ _ +0 J A ...dE, we get: 



lim / dE- / dk 



71 



■\(E 2 (k)-z)\* 



\F v (k)\ 2 dk. 



(4.33) 



{keT 2 :_B 2 (k)GA} 



So we are left with the second term of the r.h.s of Q4.31|) . For every j G {— N'J , N' q }, 
k = (ki, k 2 ) G fj, and E G A, we have: 



KmF £+ie (k) = y*(x,k;E)/(x), 



(4.34) 



where 



*(x,k,S) 



E 2 (k)-E 



„2i7rk-x 



(4.35) 



+ t m (k 2 — ma; E + iO)G^(xi; E + iO + cos27r(A;2 — ma))e 



2iir(k2—ma)x2 



m=0 



which in particular corresponds to [(E 2 (k) — £^(/c 2 ))/j(/c2)] _1 ^/ S j(x, k 2 ) for E = Ej(k 2 ). The 
limit ( 4.34|) is also uniform in k = (ki,k 2 ) G Uj and in E G A. Applying again the same 
operation: lim^o f A ...dE to the jth term of the sum in r.h.s of [4.31| , we get 



lim / dE- / dk 



71 



'\{E 3 {k 2 )-z)\< 



\FJk) 



\Fj, s (k 2 )\ 2 dk 2 . 



(4.36) 



Relations ( |4.33| ), and ( |4.36| ) imply assertions (i) and (ii) of the theorem for the case of a function 
/ with compact support. The proofs of these assertions for an arbitrary function / G Z 2 (Z 2 ), and 
the proof of assertion ii) require standard means of spectral theory (see the proof of Theorem 



Proof of Proposition [O] According to (|4.1|) , we can write equation ( |3~? ) for a = p/q as 

<? 

P q {k 2 , E) = a" ]]b{k 2 + l/q, E + iO) = 1, (4.37) 

where 



i=i 



b(k 2 ,z) 



gG^(0,z + cos27rA; 2 ) — % 



gG$\0,z + cos27rk 2 ) + i 



(4.38) 



Since the product in the l.h.s. of the equation [4.37| is periodic in k 2 with period l/q, its 
solutions are also periodic in k 2 with period l/q, and we can restrict ourselves to the interval 
[1/2 - l/(2g),l/2 + l/(2g)). By Lemmas [H| - |5j we have |6(jfe 2 ,«)| < 1, e T, z G C, 
thus the band equation ( f4.37| ) admits a solution if and only if the modulus of each factor 
b(k 2 + l/q,E + «0), I = 1, q in its l.h.s is 1. Hence, we can write the representation 



b{k 2 ,E) = exp{27ii(j)(k 2 ,E)}. 
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In what follows we will consider the case where E is positive (the arguments for negative E are 
similar and will be omitted). In this case we have from (|3.37|) : 

G 1 \o,E + tO) = -^=, E>1, 
and we can choose the phase <p(k 2 ,E) as 



(j>{k 2 ,E) = l/vr arctan ((1/g) y/(E + cos 2nk 2 ) 2 - l) . (4.39) 

For any k 2 G [1/2 — l/(2q), 1/2 + l/(2q)) is a non-negative and an increasing function of 
E > 1 — cos27r/c2, satisfying the inequalities: 

0<<j){k 2 ,E) < 1/2. 

The above formulas show that the l.h.s of equation ( [4.371 ) is rea l analytic in the domain 

{(k 2 , E):k 2 e [1/2 - l/(2g), 1/2 + \/{2q)), E > 1 - cos2vrfc 2 }, 

hence solutions of the equation, if they exist, are real analytic in k 2 (notice that here the 
condition E > 1 — cos27rA;2 is equivalent to the conditions E > 1 — cos 1n{k 2 + la), V/ = 1, ..,q). 

We will use ( [4.37 ) in the form 

<5> q (k 2 ,E) - qcu = (modi), (4.40) 

where 

q q-l 

<$> q (k 2 , E) = J2 <P( k 2 + la, E) = J2 <P( k 2 + la, E). (4.41) 

1=1 1=0 

For any fixed k 2 G [1/2 — l/(2g), 1/2 + l/(2g)), $ g (fc 2 , i?) is a positive and an increasing function 
of £7 > 1 — cos27r/c 2 , bounded by q/2. 

Fix now g and uj and denote by a u the integer part of the minimum 

min & (k 2 , 1 — cos 2nk 2 ) — quj. 

fc 2 6[l/2-l/(2g),l/2+l/(2 g )) 

For a fixed integer j denote by Ej(k 2 ) the energy such that 

$ q (k 2 , Ej(k 2 )) -quo = a UJ + ] (4.42) 

and denote by Vj the set of k 2 G T such that ( |4.42j ) is satisfied. The sets Vj, j = —N^, ...,N' 
form an increasing family of the closed subset of T. For all j larger than some jo , coincides 
with T. 

Hence Ej is the j-th energy band function and Rani^- is the j-th surface energy band. It is 
clear that the maximum value N' q of j for which such a solution exists, is such that N' < q/2. 

Since Q q (k 2 , 1 — cos2Tck 2 ) — qu> is analytic in k 2 G T, it may exist a discrete set T>'a of k 2 G T, 
for which <& q {k 2 , 1 — cos 27tk 2 ) — qu) is equal to the integer + j. Numerical experiments show 
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that for small q there are at most two values of k 2 in the interval [1/2 — l/(2q), 1/2 + l/(2g)) 
for which this event occurs, so the number of points in T>[ is 2q and the other Vj are empty. 

We have proved that if k 2 € Tj = T>j \ Vj, Ej(k 2 ) exceeds 1 — cos27rA; 2 , then we have for all 



Ej(k 2 ) > 1 — cos27r£; 2 > — cos27r/ci — cos27r/c 2 = ^(k), 



(4.43) 



i.e. the separation property (|4.8|) between the band of the volume states and the surface bands. 
Let us now discuss separation between the surface bands E_ N n, ...E N >. We will use the relation 



^ q (k 2 ,E j+1 (k 2 ))-^ q (k2,E j (k 2 )), 



(4.44) 



implied by fl£40P . 

Consider first the energy range E > e m for some e m > 2. It follows from ( |4.40|) that the 
maximum energy E q for which the equation is soluble is finite (this is the upper edge of the 
spectrum of the operator H for a given q). Hence the partial derivative 



9 



E + cos27r(fc 2 + l/q) 



BE 

satisfies the inequalities: 



9-1 

1/7F S g 2 + [ E + cos 2vr(A; 2 + //g)] 2 - 1 ^(E + cos 2tt(A; 2 + //g)) 2 - 1 



(4.45) 



< 



9$ 



5^ tt(7 ^/(e,,, - 2)e m ' 
This bound and ( [4.44 ) lead to the relations 

r E ^ k2) d%{k 2 ,E) 

iE(k 2 ) 9E 



(v. 



>\-i 



< oo. 



1 



dE < (E j+1 (k 2 ) — Ej(k 2 )){r]'] 



(4.46) 



implying the separation property (|4.9| ) in the case where Ej(k 2 ) > 2. 
In the case, where 

< 1 - cos2ttA; 2 < Ej(k 2 ) < 2, k 2 G [1/2 - l/(2g), 1/2 + l/(2g)), 



the r.h.s of ( |4.45| ) can be infinite because of the contribution of the first term (for E — 1 — 
cos27rfc 2 ), and of the second term (for E — 1 — cos27rfc 2 , and k 2 = 1/2 — l/(2g)) or of the 
(q — l)th term (for E = 1 — cos27r/c 2 , and fc 2 = 1/2 + l/(2g)). Since, however, each term in the 
phase ( [4.41|) is non-negative and 



max <t>(k 2 , 2) 

*ae[i/2-l/(2ff),i/2+i/(2ff)) 



is strictly less than 1/2, the contribution of these terms in the difference (f4.44|) is bounded from 
above by 2<fr q < 1, and we obtain from ( f4.44|) the inequality 



< 1 



< 



34 



^ d%(k 2 ,E) 



dE 



dE, 



where §> q (k 2 ,E) is the sum in ( |4.41| ), in which the terms corresponding to I = and to I = 1 
if k 2 E [1/2 - 1/(2?), 1/2), and to I = q - 1 if k 2 G [1/2, 1/2 + l/(2q)) are omitted. It is easy 
to check that the partial derivative of & q (k 2 , E) with respect to E is bounded from above by a 
constant (r/q) 1 < oo. This leads to the bound ( |4.46|) in which (77' ) 1 is replaced by (rj^)' 1 and 
$ g by Q q . The obtained bounds imply the separation property (|4.9| ) with r\ q = mm{r] q , r/^} . 
Proposition [O is proved. 



Remark. It can be seen from the proof above that the distance between the bands increases as 
\j\ increases. Besides, the distance between the two first bands is of order 0(l/q) when q is 
large. 

Denote from now on the operator of ( [L.l|) - ( |1.4|) as H a . We conclude this section by discussing 
correspondence between the spectrums of the operators H a with an irrational number a and 
with its rational approximations p n /q n '- 

lim V — = a. (4.47) 



It is easy to prove, by using the basic formula ( |2.16| ) for the resolvent, that under condition ( 4.47| ) 



H Pn / qn converges to H a in the strong resolvent sense. Hence, according to general principles [[L8|], 
the spectrum a(H Pn / qn ) is upper semi-continuous in n in the limit ([4.471) . Here is a statement, 



that gives a more detailed behavior of a(H p / q ) for large q. Recall, that a(H p / q ) is the union 
of the interval [—2, 2] and of N' q ' + N' surface bands, part of which can belong or intersect the 
interval [—2, 2]. 

Theorem 4.7. Assume that q sufficiently large. Then there exists at most one negative surface 
energy band above E = —2 and at most one positive surface energy band below E = 2. These 
bands, if they exist, have the width of order 0(1/ q 2 ) as q — > 00. The width of the surface energy 
bands lying in (—00, —2) and in (2, +oo) ; are of order 0(exp{— const • q}) as q — > 00. 

Proof. We start from the dispersion equation fl4.41|) - ( [4.40| ) for the surface energy bands. Since 
the function Q q (k2,E) has period 1/q in k 2 , its Fourier series is: 

%(k 2 ,E) = J2%AE)e- 2mk ^ n , 

where 

%, n (E) = q dk 2 e 2mk ^ n <P{k 2 + l/q, E) = q f dk 2 e 2 ^ n (f)(E , k 2 ) := qj> qn (E), 
Jo Jo 

and (j) qn is the gn-th Fourier coefficient of the function </>(., E). Hence 

%(k 2 ,E) = q(ME) + UE)e- 2 ^ + ...). (4.48) 

According to ( |4.39|) , the function <f)(-,E) is analytic for \E\ > 2, thus its Fourier coefficient 
4> q (E) is of order exp{ — const - q} as q — » 00. In addition, formula ( 4.40|) implies the relation 



dEj _d$ q 



dko dko \ BE 
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It follows now from Q4.48| ) and from the exponential decay of the Fourier coefficient & q ^ n (E) 



that the upper bound for the derivative — — is of order O(q 2 e~ constq ), while the lower bound 

ok 2 

for — — — , which is reached for the highest energy band E^i(k 2 ), is of order 0(q). Thus the 
oE q 
dE- 

derivative — — is of order 0(ge _const ' 9 ). Since Ej(k 2 ) is periodic in k 2 with period 1/q, then 

dk 2 

denoting respectively by E 1 J hax , and by E™ m the maximum and the minimum of the j-th band 
function Ej, we see that \E™ ax — EJ l%n \ is of the order exp{— const • q] if EJ 1111 > 2. 

Let us fix k 2 G [1/2 — l/(2q), 1/2 + l/(2q)). To see how many bands are in between the lowest 
possible energy E = 1 — cos2nk 2 and the energy, E = 2, let us calculate 5$ = & q (k 2 ,2) — 
& q (k 2 , 1 — cos27r/;; 2 ). We have: 

5* = 1/tt g f g- — , E + cosMfa + Z/g) 

^Ji- M k 2 9 2 + (E + cos2tt(A; 2 + Z/g)) 2 - 1 yj (E + cos27r(A; 2 + l/q)) 2 - 1 

Performing the integration for the different values of I and summing respective contributions 
we obtain that 5$ is of the order (1/q) logq as q — > 00, thus 5$ — >• as g — ► 00. Remembering 
that for each fc 2 the energy of a band corresponds to an entire value of $ g + qu, we deduce that 
for large q there is at most one band in the interval to the left of 2. Since the minimum of E\ 
for k 2 6 [1/2 — l/(2q), 1/2 + l/(2g)) is larger than 1 + cos7r/g, the width of any band, lying 
inside the interval [1 + cosn/q, 2), is bounded by n 2 /2q 2 . 

It can also occur that E™ ax > 2. In this case, the same argument as above show that the part 
of the energy band in (2, 00) is exponentially small in q. Thus the total width in that case is 
at most of the order 1/q 2 . 

Remark. The assertion of the theorem can be interpreted as a kind of continuity of the spectrum 



with respect to the limiting transition ( 4.47 ). Indeed, according to the theorem, the width of the 
surface bands of H Pn / qn , lying outside the interval [—d, d], is exponentially small in q n — » 00. It 
can also be shown that the gaps between these bands are of the order l/q n . This is in agreement 
with the "limiting" form of this part of the spectrum of H a for irrational a's, satisfying the 
Diophantine condition (|1.6|) . Indeed, according to [ I§f , the spectrum of H a in this case is pure 



point and dense on R \ [—d, d\. Here is one more manifestation of this continuity. 

Recall that according to |L9] the eigenvalues of H a outside [— d, d] are indexed by x 2 G Z d2 , and 
for each x 2 G Z d2 the eigenvalue E x , 2 is the unique solution of the equation 

f(E X2 ) = a-x 2 + u (modi), (4.49) 

where / : R \ [—d, d] — > R is the monotone increasing function, defined for E > d as 

f(E) = -- f rfA; 2 arctanf 91o (k 2 , E)Y\ (4.50) 

7T Jjd 2 \ J 

or, in view of ( |2.23| ) and ( p.37| ), and for d 2 = 1 

f(E) = -- dk 2 arctan g' 1 yftjt + cos2ttA; 2 ) 2 - 1. (4.51) 

7T Jjl 
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On the other hand, we can write the band equation ( J4.40 ) as 



— $ q {k 2 ,E) = — +cu 



(4.52) 



for some integer I. Choosing I in the form I = p n X2 + mq n for some integer m, we can write the 
last equation as 



— $ qn (k 2 ,E) = —x 2 + uj. 

q n q n 



(4.53) 



Recalling now the expression (|4.41 ) for the function Q q {k,2, E), we conclude that for the limiting 
transition fl4.47|) and E > 2 the equation ( |4.40|) , defining the surface bands of H Pn / qn outside 
[— d, d], converges to the equation (|4.49| ), defining the all eigenvalues of H a for a Diophantine 
a outside [— d, d). 



5 Auxiliary Facts 



We present here useful facts on the Green function ( [2.18 ) of the //-dimensional Laplacian and 
on related quantities. 



Lemma 5.1. Let G^ipc — y;z), x, y £ 7j u ,^sz ^ 0, be the Green function (\2.1S ) of the v- 
dimensional Laplacian Write 

G%\0; z) = R u {z) + il v (z), R u , l v e R. (5.1) 

Then 

(i) for any e > 0, and Eel 

\R V (E + ie) | < oo, <I u (E + ie)< oo; (5.2) 

(ii) the limits R u (E + iO) and I v (E + iQi) exist for \E\ ^ v , satisfy inequality ( \5. Sj ) for \E\ < u, 
and 1 U {E + iO) = if and only if \E\ > v. 



Proof. The part (i) of the lemma follows from the integral representation (2.5). It is also easy 
to prove that the limits R V {E + zO) and l v (E + iO) exist and are finite for \E\ ^ v (in fact, 
for v > 3 they are finite even for \E\ = v, see Lemma |5.4j below). Thus we have to prove that 



l v [E + iQ) is strictly positive for \E\ <v. By using fl2.18|) , it easy to show that for v = 1 
Tr^IxiE + iO) := P i{E) = | 



[1-E 2 )- 1 / 2 , \E\ < 1, 
0, \E\ > 1, 



and that tt 1 I u {E + iO) is the uth convolution of p\. These two observations imply that 
ix~ x l v {E + iO) is strictly positive if \E\ < v, and is zero for \E\ > v. Lemma is proved. 
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Lemma 5.2. Let 70 (z) be the operator in l(Z d2 ), defined as 

7o(X 

and 

91o(z) - i 



Pz^Gq (z)Pjd 2 , d 2 < d, 



b(z) 



97o{z) + i 



(5.3) 



Then the operator 70 (z) + i is invertible for $sz > ; and the operator b(z) is a contraction for 
Sz ^ 0: 

||6(z)||<l. 

Proof. According to (|2.18|) and Q2.22| ) 70 (z) is the a convolution operator in Z 2 (Z d2 ) and its 
symbol 70(^2; z) satisfies the inequality: 37o(/c2; z) > 0, ^sz > 0. Since the symbol of 70(2) + i 
is 70(^2; z) + i, we have that (70(^2; z) +i\ > ^(70(^2; 2) +i) > 1- Hence 70(2) + i is invertible 
and ||(7o(2) + < 1. 

The operator b(z) is a rational function of Jq(z), thus it norm can be found as 

||6(z)|| = sup \b(k 2 ;z)\. 

k 2 eT d 2 



By using (|5.1|), we obtain that 

b(k 2 ,z) 



K + (I* + I) 2 



(5.4) 



where i?,, and are defined in (15. II) . This formula and Lemma 15. II lead to (157 



Lemma 5.3. Let b(k 2 ; z) be defined by $i.23j . Then 
(1) \b{k 2 ; E + z0)| < 1, V£ G E; 

(mJ /or any 7 > and < a 1 — 7 taere exists an open set K^(E) C T d2 , such that 

b(k 2 ;E + iO) < 1, fc 2 G K 7 (E). 



(5.5) 



Proof. The part (i) of the lemma follows from Lemma |5.1| , and from ( |5.4j ). To prove assertion 
(ii) we have to find that for any 7 > and \E\ < d — 7 there exists an open set K y (E) such that 
for k 2 G K^(E), \E — E d2 {k 2 )\ < d\. Then I dl (E + iO) will be strictly positive and b(k 2 ;E + iO) 
will be strictly less then 1 in view of ( |5.4j ). Since E d2 is a continuous function in k 2 on T d2 , 
varying between — d 2 and d 2 , respective open set K y (E) always exists if \E\ < d. Lemma is 
proved. 

Lemma 5.4. Let Gq (x.;z) be the Green function of the is -dimensional Laplacian and g > 0. 
Then the expression 



G { o\x;E + iO) 
gG ( o\0;E + iO)+i 



(5.6) 



is bounded in x G 27 and in E G 
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Proof. Consider first the one- dimensional case v = 1. Then it follows from ( [4.14 ) that the 
expression ( |5.6|) is 

gj??(E+iO)|2;| 



g + sinr](E + iO) ' 

and, according to ( |3.38| ) - ( |3.39| ), the modulus of the last expression is bounded by g~ x . 
For v > 2 we will use the integral representation of Gq , (x; z) of (|2.18| ), valid for ^sz > 0: 



dte 4 ** Y[j Xl {t)e 



inxi/2 



(5.7) 



where x = {^}f =1 , and J n (i) is the Bessel function of the order n: 

-2tt 



J n (t) 



1 

2^ 



ini)+it sin 



The representation follows easily from (|2.5| ), and from the identity 

/■oo 

(A - z) _1 = i I dte- it(x - z \ A 6 1, 82; > 0. 
Jo 

By using the asymptotic formula 



(5. 



Jn(t) 



we find that v > 3 Gq (x; + iO) is bounded in x and in £7. Since, in addition, |gG yJ (0; + 
i0) + i| > tfQG^O; £ + zO) + 1 > 1 (recall that in view of (gig) SG^O; 2) is nonnegative), 
we obtain the assertion of the lemma for z/ > 3. 



2 V /2 , 

cos,f 



00, 



(5.9) 



Thus, we are left with the case z/ = 2. By using again (|5.7| ) and (p.9|), we find that Gq^(o;; + 
iO) is bounded in x and in E everywhere except \E\ = 2, and that in a sufficiently small 
neighborhood of E = 2 



G^(x; £ + iO) = A{x) log \E - 2| + £±(x) + 0(|£ - z|), £ 



0. 



where A(x) and 5(x) are bounded in x, A(0) 7^ 0, and B±{x) correspond to sign(_E — 2). The 
same asymptotic representation is valid in a neighborhood of E = —2. This shows that the 
ratio G^\x; E + iO) / G^\0; E + iO) is bounded and continuous in E G K for any x G If . In 
addition we have: 



G [ u\0;E + iO) 



gG i l/ \0;E + iO)+i 



1 



g + i 



< 



G ( o\0;E + iO) 
1 



9 ~ 3 



-1 
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because 

-l 



-3 



G { o\0;E + iO) = %G { l/> (0-E + iO)/\G { Q l/> (0;E + iO)\' 2 > 0. 



Lemma is proved. 

Lemma 5.5. The expression 

g^oih, E + iO) +i 
is bounded in E 6R, k 2 G T d2 , and x G Z dl . 

Proof. According to fl2Jj), 7(^3, z) = Gj dl) (0; z - E d2 (k 2 )). Hence, we can apply Lemma |5.4| 
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